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19.  ABSTRACT  (continued).... 


loading  chains  and  islands  of  quasistatic  subassemblies,  and  (4)  the  relationship  between  the  number  of 
intergranular  contacts  and  macroscopic  strain. 

The  behavior  of  the  frictional  material  is  shown  to  be  very  complex,  involving  a  number  of 
different  physical  processes.  To  gain  a  better  understanding  of  the  micro  mechanical  behavior  a  (VHS 
compatible)  video  tape  was  generated  of  the  response  of  the  particle  assembly  throughout  the  numerical 
simulation.  A  review  of  the  response  of  the  frictional  material  during  the  test  shows  the  development 
of  inhomogenities  associated  with  the  acretion  of  particles  to  quasistatic  islands  of  particles  in  the  sam¬ 
ple. 


The  experimental  portion  of  the  research,  described  in  Section  III  of  this  report,  has  been 
directed  toward  developing  a  high  pressure; high  precision  test  facility  capable  of  exercising  a 
moderately  large  cylindrical  specimen  (3.0  in.  diameter  by  7.0  in.  height)  to  stresses  of  the  order  of 
1.0  kbar  (15  ksi)  in  a  conventional  triaxial  test.  The  objective  of  this  effort  was  to  build  a  device  with 
sufficient  sensitivity  to  permit  an  experimental  examination  of  inelastic  material  behavior  over  a  wide 
range  of  conditions  where  a  number  of  important  theoretical  stability/instability  issues  have  been 
raised. 


While  it  was  not  possible  to  construct  the  facility,  conduct  a  comprehensive  test  program,  and 
coordinate  the  test  results  with  the  numerical  analyses  within  the  two-year  grant  period,  the  facility  has 
been  completed  and  tested  at  a  range  of  0.1  kbar  where  results  can  be  compared  with  existing  test  data. 
Based  on  such  comparisons,  the  test  apparatus  appears  to  be  quite  sensitive  and  able  to  track  specified 
load  or  strain  inputs  at  the  level  of  approximately  0.1%  of  maximum  values.  The  machine  is  designed 
to  investigate  quasistatic  stress-strain  behavior;  however,  the  machine  is  driven  by  a  high  speed  micro¬ 
computer  driven  control  system  and  is  limited  only  by  the  response  times  of  the  stepper  motors  which 
drive  the  hydraulic  loading  system.  It  is  believed  that  the  experimental  test  facility  is  capable  of 
exercising  moderately  large  samples  of  a  frictional  material  over  a  wide  range  of  triaxial  loading  condi¬ 
tions  where  theoretical  issues  regarding  material  stability  can  be  explored  experimentally. 

The  results  of  the  experimental  study  will  be  forthcoming  as  addendums  to  the  final  report  as 
they  become  available. 
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I.  INTRODUCTION 


The  research  presented  in  this  report  is  the  result  of  a  two-year  "Combined  Experimen¬ 
tal  and  Numerical  Investigation  of  Inelastic  Behavior  of  Frictional  Geologic  Material,"  sup¬ 
ported  by  AFOSR  Grant  #86-0290.  "^The  objective  of  this  research  was  to  develop  a  coordi¬ 
nated  experimental/analytical  capability  for  investigating  the  basic  physics  of  granular  geolo¬ 
gic  materials  emphasizing  their  microphysical  behavior.  The  research  is  motivated  by  the  fact 
that  the  inelastic  behavior  of  geologic  materials,  except  for  very  simple  and  thoroughly  exam¬ 
ined  loading  conditions,  has  been  difficult  to  predict  using  classical  continuum  mechanical 
theories.  A  major  conflict  between  theoretical  and  experimental  results  arises  where  labora¬ 
tory  tests  demonstrate  material  behavior  which  cannot  be  captured  using  existing  theoretical 
models  for  geologic  materials.  The  primary  issues  have  been  centered  around  nonassociated 
flow,  where  inelastic  flow  is  observed  to  differ  significantly  from  theoretically  accepted  asso¬ 
ciated  flow  (perpendicular  to  material  yield  surfaces.)  According  to  conventional  continuum 
mechanical  theory  such  materials  will  be  unstable  when  exercised  along  certain  strain  or  stress 
paths.  Yet  in  laboratory  tests  the  materials  may  or  may  not  be  stable. 

Recognizing  that  continuum  mechanical  '^stress-strain*  models  are  fundamentally 
phenomenological  in  nature  the  current  research  is  directed  toward  investigating  geologic 
material  behavior  at  the  micromechanical  level  using  both  numerical  and  experimental 
methods.  JUie  results  of  a  numerical  analysis  of  micromechanical  behavior  of  granular  materi¬ 
als  are  contained  in  Section  II  of  this  report,  which  is  a  condensed  version  of  a  Ph.D  disserta¬ 
tion  of  one  of  the  co-authors  (Issa).'^ection  III  is  a  description  of  a  high  pressure,  high  preci¬ 
sion  test  facility  constructed  during  the  contract  period.  The  facility  is  intended  to  be  used  for 
high  precision  triaxial  tests  of  granular  materials  over  a  wide  range  of  stress  levels  where 
many  issues  regarding  nonassociated  flow  and  the  related  issued  regarding  stability  of  geolo¬ 
gic  materials  can  be  resolved  experimentally.  \ 
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As  will  be  seen  both  the  numerical  and  experimental  programs  are  still  under  way  and 
the  two  activities  are  not  completely  merged  at  present.  The  numerical  analysis  provides  very 
interesting  and  realistic  simulations  but  has  been  restricted  to  small  assemblies  of  (200-300) 
grains  due  to  computer  budget  limitations.  This  restriction  has  also  prevented  a  full  examina¬ 
tion  of  a  number  of  questions  regarding  parametric  sensitivities  (geometry  of  grains,  contact 
friction,  grain  elasticity  and  fracture  strength,  etc.)  The  experimental  test  facility  required  a 
longer  time  to  complete  than  initially  estimated.  However,  the  facility  is  now  completed  and 
high  pressure  testing  is  in  progress.  The  device  appears  to  have  met  all  the  high  pressure/high 
sensitivity  requirements  set  at  the  beginning  of  the  project.  Given  these  new  capabilities  for 
conducting  numerical  and  experimental  analyses  of  frictional  granular  materials,  it  is  expected 
a  number  of  new  findings  regarding  the  fundamental  behavior  of  geologic  materials  will  be 
produced  as  research  continues.  These  results  will  appear  as  Addendums  to  this  Final  Report. 
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n.  NUMERICAL  ANALYSIS  OF  MICROMECHANICAL  BEHAVIOR 
OF  GRANULAR  MATERIALS 


H.l.  Introduction 

Perhaps  the  most  striking  characteristic  of  granular  matter  is  its  dual  nature  which  lies 
between  a  disjoint,  discrete  material  and  a  continuum.  For,  although  the  individual  particles 
of  the  medium  are  solid,  the  particles  are  only  partly  connected.  Thus  the  basic  structure  of 
the  material  varies  greatly  with  loading  and  material  behavior  can  range  from  solid  to  fluid 
with  or  without  localized  secondary  physical  processes,  such  as  fracture,  and  shear  slip  band 
formation. 

A  granular  medium  may  be  defined  as  an  aggregate  of  discrete,  solid  grains  which 
occupy  a  nominal  volume.  The  grains  are  free  to  displace  with  respect  to  their  neighbors  and 
the  space  between  the  particles  may  contain  fluid,  interstitial  particles  or  be  empty.  The 
mechanical  behavior  of  an  assembly  of  granular  materials  is  a  complex  integrated  effect 
defined  by  the  geometry  of  the  structure  and  physics  of  contact,  and  in  particular  on  the 
number  and  strength  of  the  contact  bonds  which  are  themselves  a  consequence  of  the  size, 
shape,  roughness  and  strength  of  the  particles.  All  these  factors  are  incorporated  into  the  con¬ 
cept  of  structure.  The  constraint  on  the  motion  of  the  individual  particles  can  be  either  negli¬ 
gibly  small  or  complete  i.e.,  ranging  from  a  suspension  of  particles  to  multicontacted  particles 
forming  a  material  with  continuous  force-deformation  characteristics,  i.e.,  a  continuum. 
Perhaps  the  most  significant  feature  of  granular  materials  is  that  inelastic  deformation  is 
brought  about  by  mutual  sliding  of  the  particles  (intergranular  deformation),  in  contrast  to  the 
deformation  of  the  individual  panicles  (intragranular  deformation)  or  dislocation  transport  typ¬ 
ical  of  a  continuous  medium.  The  ability  of  the  body  to  maintain  its  shape  to  continue  in  the 
anisotropic  state  of  stress,  is  a  consequence  of  internal  friction. 


The  differences  in  the  behavior  of  granular  materials  when  compared  with  the  behavior 
of  continuous  materials  are  large  enough  to  consider  the  mechanics  of  granular  materials  to  be 
an  autonomous  branch  of  mechanics  which  is  not  directly  related  to  the  mechanics  of  solid 
continuous  materials.  Granular  materials  exist  in  a  wide  variety  of  technological  fields.  The 
oldest  is  soil  mechanics  with  its  problems  of  foundations,  retaining  walls,  hill  side  stability 
and  excavation.  Such  problems  have  been  treated  on  a  more  or  less  scientific  level  for  at 
least  two  centuries.  The  dynamical  response  of  granular  material  has  been  investigated  in 
connection  with  the  propagation  of  seismic  waves,  petroleum  prospecting,  and  in  the  past  25 
years,  with  response  of  geological  materials  to  blast  and  ground  shock.  The  internal 
configuration  of  granular  matter  is  found  to  be  important  in  a  variety  of  fields  such  as  in  the 
manufacture  of  ceramics,  cement  and  mortar  products,  grinding  wheels,  powder  metallurgy 
forming  processes  and  composite  materials. 

The  behavior  of  granular  matter  has  attracted  the  attention  of  researchers  since  the  19th 
century,  yet  there  still  remain  many  unanswered  questions  especially  concerning  inelastic 
deformation  and  flow.  Many  theories  have  been  developed  on  the  macroscopic  level  by  treat¬ 
ing  the  granular  assembly  as  a  continuum.  Today  large  finite  element  and  finite  difference 
computer  programs  use  such  material  representations  or  models  to  simulate  the  behavior  of 
granular  matter.  The  credibility  of  these  simulations  resides  directly  on  the  complex 
mathematical  models  of  granular  materials  designed  to  predict  material  behavior  in  realistic 
(arbitrary)  loading  conditions. 

Nearly  all  of  the  mathematical  models  of  frictional  materials  are  based  on  nonlinear 
mathematical  continuum  mechanics  theories  ranging  from  simple  curve  following  incremental 
elasticity  models  to  complex  hybrid  incremental  hypoelasticity  nonassociated  flow  plasticity 
models  [1,2, 3, 4, 5, 7].  Current  state  of  the  art  continuum  models  for  frictional  materials  are 
phenomenological  in  nature,  that  is,  constructed  to  reproduce  macroscopic  behavior  observed 
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in  experiments,  with  little  attention  paid  in  their  development  to  the  physics  of  the  micros¬ 
copic  behavior  of  granular  interaction. 

This  approach  has  been  successful  in  predicting  the  behavior  of  metals  where  contin¬ 
uum  models  are  based  on  incremental  plasticity  theory.  It  is  not  surprising  therefore  that  the 
historical  development  of  continuum  models  for  granular  materials  followed  the  same 
approach,  nor  that  as  a  first  step  the  successful  metal  plasticity  models  were  extended  to 
account  for  pressure  sensitive  frictional  behavior. 

Unfortunately  the  extension  of  metal  plasticity  based  models  to  granular  media  has  led 
to  difficult  problems  which  have  yet  to  be  solved.  Generally,  phenomenologically  based  con¬ 
tinuum  models  are  able  to  reproduce  most  test  data  reasonably  well,  especially  for  tests  where 
the  material  is  subjected  to  (nearly)  monotonically  increasing  loads  to  the  point  of  failure. 
For  tests  involving  cyclic  loadings,  even  the  best  of  today’s  models  have  difficulty  in  captur¬ 
ing  qualitative  effects  reliably  or  in  loading  regions  different  than  those  used  in  the  model 
development.  In  other  words,  the  capability  of  current  models  to  predict  behavior  of  frictional 
materials  to  arbitrary  loads  is  a  major  problem. 

To  make  matters  worse,  the  development  of  modem  continuum  models  capable  of 
accurately  simulating  test  data  has  introduced  a  number  of  perplexing  issues  of  a  theoretical 
nature  [1].  Certainly,  any  mathematical  model  should  not  violate  the  laws  of  thermodynamics 
(for  example,  by  generating  energy  in  closed  loadings  cycles)  nor  should  the  model  predict 
discontinuous  behavior,  where  two  arbitrarily  close  finite  stress  paths  produce  finite 
differences  in  strains,  unless  such  a  process  occurs  physically  as  by  say  a  material  failure.  A 
third  requirement  is  that  the  mathematical  relationships  should  be  unique  in  the  sense  of 
predicting  a  one  to  one  correspondence  between  incremental  changes  in  stress  and  strain. 
Finally,  the  model  should  not  indicate  unstable  material  behavior  if  in  fact  the  behavior  of 
physical  material  is  stable,  and  vice  versa. 


In  an  effort  to  account  for  realistic  material  behavior,  complex,  but  theoretically  "safe" 
models  ba„cd  on  associated  flow  incremental  plasticity  have  been  developed  [4,5,6].  These 
material  models  can  be  made  to  obey  all  the  basic  thermodynamic  and  continuity  requirements 
imposed.  Unfortunately,  such  models  do  not  give  a  realistic  representation  of  ordinary  materi¬ 
als  because  inelastic  strain  increments  developed  when  the  material  is  exercised  in  shear  are 
not  perpendicular  to  the  yield  surface,  a  characteristic  phenomenon  known  as  nonassociated 
plastic  flow  for  frictional  granular  materials  [35]. 

A  different  attempt  to  treat  the  observed  behavior  of  granular  materials  has  evolved 
from  high  pressure  physics,  where  the  basic  information  describing  material  behavior  is  in  the 
form  of  a  pressure/volume  relationship  [1].  For  problems  in  the  low  pressure  region  shear 
strength  becomes  important.  To  accommodate  this  effect,  shear  strength  was  appended  to  the 
pressure/volume  relationship  using  Prandtl/Reuss  (volume  constant)  plasticity  relations  to 
describe  shear  induced  inelastic  behavior.  These  models  are  hybrid  in  the  sense  that  the 
pressure/volume  effect  is  partially  decoupled  with  the  pressure  dependent  shear  behavior. 
These  material  models  can  be  made  to  fit  test  data  significantly  better  than  the  associated  flow 
models.  Unfortunately  these  models  can  be  shown  to  violate  some  of  the  theoretical  require¬ 
ments  from  basic  physics,  namely  the  requirement  against  energy  generation  and  sometimes 
the  continuity  requirement  as  well. 

In  summary,  it  can  be  said  that  while  associated  flow,  hybrid  elasto-plastic  and  nonas¬ 
sociated  flow  plasticity  type  models  may  give  reasonable  fits  to  test  data  and  adhere  reason¬ 
ably  well  to  basic  principles,  no  single  model  is  able  to  capture  all  observed  behavior  of 
granular  materials  and  at  the  same  time  meet  all  theoretical  requirements.  This  is  in  spite  of 
more  than  thirty  years  of  research  on  the  problem  of  developing  a  satisfactory  mathematical 
representation  for  granular  materials. 


6 


Ther  .ore,  this  research  is  directed  towards  attempting  to  better  understand  the  basic 
physical  processes  which  are  dominating  the  behavior  of  granular  matter.  These  processes  are 
all  believed  to  be  at  a  localized  interparticle  level.  That  is,  macroscale  "stress-strain"  material 
behavior  is  a  consequence  of  the  physics  of  contact  and  friction  between  particles  of  granular 
materials.  If  this  behavior  can  be  clearly  understood  then  it  is  believed  that  it  will  be  much 
easier  to  establish  a  theoretically  sound  and  accurate  continuum  model. 

Most  of  the  early  analytical  work,  in  an  attempt  to  refine  the  traditional  continuum 
models,  took  natural  step  in  the  direction  of  replacing  a  continuous  medium  by  one  in  which 
account  is  taken  of  forces  transmitted  between  individual  particles.  Initial  work  along  these 
lines  approximated  the  particles  by  rigid  spheres. 

Hara  [8]  seeking  to  predict  the  behavior  of  the  carbon  microphone,  studied  the  propa¬ 
gation  of  compressional  waves  through  an  aggregate  of  like  elastic  spheres  arranged  in  simple 
cubic  and  face-centered  cubic  arrays,  replacing  the  granules  by  a  series  of  mass-spring  sys¬ 
tems  in  series,  and  computing  the  stiffness  of  each  spring  on  the  basis  of  Hertz’s  theory  of 
contact. 

Proceeding  in  a  similar  vein,  Gassman  [9]  approximated  the  earth’s  crust  by  an  array 
of  like  spheres  in  a  hexagonal  close  packing.  When  subjected  only  to  its  own  weight,  and 
when  the  stiffness  at  the  individual  contacts  was  computed  on  the  basis  of  Hertz  theory,  this 
array  was  shown  to  have  elastic  constants  appropriate  for  a  continuum  with  transverse  iso¬ 
tropy.  The  abandonment  of  the  regular  lattice  served  as  the  starting  point  for  a  study  by 
Brandt  [10]  who  considered  a  model  consisting  of  several  sets  of  like  spheres,  the  spheres  in  a 
given  set  fitting  in  the  interstices  of  the  preceding  set. 

The  inclusion  of  tangential  contact  forces  by  Mindlin  and  Deresiewicz  [12]  in  the  solu¬ 
tions  of  elasticity  equations  produced  significant  corrections,  since  the  stiffness  of  a  contact  in 
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a  tangential  direction  is  of  the  same  order  of  magnitude  as  its  stiffness  in  the  direction  of  its 
normal.  Duffy  and  Mindlin  [11,12]  also  computed  the  incremental  contact  forces  and  incre¬ 
mental  stress-strain  relations  for  the  face  centered  cubic  array  under  initial  isotropic  and 
uniaxial  compressive  stress. 

Using  similar  static  analyses,  Verruijt  and  De  Josselin  De  Jong  [13],  Oda  and  Konishi 
[14]  have  revealed  the  existence  of  chains  of  stresses  in  the  material,  but,  due  to  experimental 
difficulties,  they  were  unable  to  bring  out  clearly  the  evolution  of  the  contact  forces  under 
overall  stresses. 

Kinematical  approaches  for  analyzing  the  displacements  of  the  particles  were  based  on 
the  study  of  the  geometrical  evolution  of  a  two-dimensional  granular  structure.  Studies  by 
Biarez  [15],  Wiendieck  [16],  Matsuoka  [17],  Oda,  et  al.  [18],  and  Chapuis  [19]  reveal  the 
existence  of  a  structural  anisotropy  which  evolves  with  deformation  and  which  is  a  kinemati¬ 
cal  demonstration  of  the  memory  of  the  material.  Another  recent  interesting  approach  by 
Christofferson  and  Nemat  Nasser  [18],  and  Oda  [18],  proposes  to  define  stresses  in  terms  of 
the  fabric  and  other  microstructural  characteristics  of  granular  mass.  Cambou  [21]  proposed  a 
description  of  the  state  of  a  granular  material  based  on  statical  analysis.  He  constructed  a 
numerical  model  to  show  that  this  kind  of  description  is  able  to  take  into  account  the  loading 
history  of  the  material.  He  concluded  that  the  so-called  hardening  of  granular  material  is  nei¬ 
ther  isotropic  nor  kinematic.  Matsuoka  [22]  proposed  a  stress-strain  model  for  granular 
materials  by  consideration  of  the  mechanism  of  fabric  change  under  shear.  He  derived  a  rela¬ 
tionship  between  strains  and  fabric  change  of  granular  materials  by  considering  two  particles 
and  the  mechanism  of  disappearance  and  generation  of  interparticle  contacts  with  a  third  parti¬ 
cle. 
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On  an  entirely  different  line  of  approach,  the  mechanics  of  granular  materials  can  be 
investigated  using  stochastic  techniques  and  modeling.  Attempts  in  this  direction  were  made 
by  Massal  [23],  Kitamura  [24]  and  Chikwendu  [25].  Recently  progress  was  made  by  explor¬ 
ing  the  kinetic  theory  for  rapid  flow  of  rough  inelastic  particles.  Ogawa  [26],  and  Savage  [27] 
have  been  partially  successful  in  applying  such  formulation  to  the  flow  of  granular  materials. 

Serrano  and  Rodriquez-Ortiz  [28]  developed  a  numerical  model  for  assemblies  of  discs 
and  spheres.  Contact  forces  and  displacements  are  calculated  for  equilibrium  conditions 
assuming  that  increments  of  contact  forces  are  determined  by  incremental  displacements  of  the 
centers  of  the  particles.  A  major  draw-back  of  the  method  used  to  solve  the  equations  is  that 
only  a  relatively  small  number  of  particles  can  be  processed. 

A  separate  type  of  approach  is  based  on  the  Distinct  Element  Method  (DEM),  a 
numerical  method  that  is  capable  of  handling  disconnected  particles  of  any  shape  developed 
by  Cundall  [29]  for  the  analysis  of  rock  mechanics  problems.  In  the  distinct  element  method, 
the  interaction  of  the  particles  is  viewed  as  a  transient  problem  with  states  of  equilibrium 
developing  whenever  the  internal  forces  balance.  Cundall  and  Strack  [30]  used  a  version  of 
discrete  element  method  named  Ball  to  investigate  two-dimensional  behavior  of  arbitrarily 
sized  circular  grains  subjected  to  a  variety  of  strain  paths. 

Bazant  [31]  recently  used  a  simplified  numerical  method  to  study  the  microstructure 
and  crack-growth  in  geomaterials.  This  method  is  somewhat  related  to  the  discrete  element 
method  as  the  frictional  interaction  is  replaced  by  a  force-displacement  relation  with  a  tensile 
strength  limit.  However  block  motion  and  the  on-off  nature  of  contacts  were  incompletely 
treated. 


The  discrete  element  method  serves  as  a  starting  point  to  developing  an  understanding 
of  the  primary  physical  processes  which  govern  the  behavior  of  granular  materials.  As  with 
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earlier  research  using  DEM  the  level  of  modeling  will  be  two-dimensional  and  at  the  micro¬ 
physical  level  where  assemblies  of  arbitrarily  shaped  and  oriented  material  particles  will  be 
analysed  to  simulate  the  behavior  of  a  real  granular  material.  For  this  purpose  a  new  program 
MASOM  (Micromechanical  Analysis  and  Simulation  of  Materials)  has  been  written  during  the 
current  research  program.  This  program  is  based  on  an  earlier  program  for  the  response  of 
block  systems  named  PROBS  [32,33,34,35],  which  in  turn  evolved  from  Cundall’s  early  RBM 
(Rigid  Block  Mechanics)  program  [29].  Although  based  on  PROBS  the  program  has  a  re¬ 
worked  solution  strategy  and  completely  new  physics  of  intraelement  and  interelement 
behavior.  In  its  current  form,  the  computer  program  is  capable  of  analysing  frictional  contact 
forces  damping  and  granular  deformation,  and  includes  an  ability  to  represent  grain  crushing 
at  high  pressure.  The  program  has  well  been  developed  to  optimize  computational  efficiency. 

The  program  is  used  to  investigate  microphysical  behavior  of  geologic  materials  by 
simulating  conventional  laboratory  tests  such  as  uniaxial  strain,  triaxial  compression  and  sim¬ 
ple  shear.  A  number  of  interesting  results  were  obtained  as  will  be  seen  in  the  following.  One 
development  which  proved  very  helpful  in  understanding  micro-mechanical  behavior  is  that  of 
a  graphics  software  capability  for  creating  a  metafile  consisting  of  video  frames  for  animation, 
i.e.,  producing  movies  of  the  test  simulation.  This  simulation  revealed  several  processes 
which  cannot  be  shown  with  conventional  time  history  plotting  or  graphical  displays  of  the 
computer  models  at  specified  times. 

II.2  Technical  Approach 

In  this  section,  several  major  features  of  a  computer  code  named  MASOM  (Micro 
Analysis  Simulation  of  Materials)  developed  during  the  course  of  the  research  are  described. 
The  program  provides  a  foundation  for  investigating  the  microphysics  of  geological  materials. 
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As  noted  above  MASOM  is  based  on  PROBS  [32-35]  and  follows  the  same  general 
methodology  for  analysing  discrete  element  systems.  As  with  PROBS,  the  MASOM  program 
tracks  the  behavior  of  an  assembly  of  arbitrarily  shaped  two  dimensional  polygons  by  means 
of  an  explicit  time  integration  scheme.  There  is  no  restriction  other  than  computational 
expense  (and  computer  memory)  on  the  size  of  the  problem  which  can  be  analysed. 

As  with  PROBS  a  number  of  different  groups  of  elements  can  be  defined,  each  with 
its  own  intraelement  and  interelement  contact  relationships  with  other  element  groups.  How¬ 
ever,  the  PROBS  computer  program,  which  was  used  to  investigate  high  velocity  impact  prob¬ 
lems,  contains  a  number  of  approximations  which  made  it  impossible  to  use  the  program  to 
investigate  low  velocity-to-quasistatic  processes  and,  in  particular,  the  basic  micromechanical 
behavior  of  granular  materials  subjected  to  quasistatic  loads  typical  of  most  conventional 
laboratory  tests. 

The  new  features  in  MASOM  which  distinguish  this  program  from  PROBS,  and  which 
make  it  possible  to  analyse  microphysical  behavior  of  frictional,  granular  materials  are 
described  next. 

II.2.1  Development  of  Consistent  Internal  Stress  Strain  Behavior  for  the  Particles 

A  consistent  methodology  for  updating  the  strain  rate  in  each  particle  is  a  basic 
requirement  for  representing  the  response  of  individual  grains.  In  PROBS  this  was  accom¬ 
plished  using  mathematical  techniques  which  could  not  be  justified  on  a  rational  physical 
basis.  Thus,  a  number  of  major  changes  were  requested.  In  the  current  procedure  three 
independent  constant  strain  rate  modes  exist  in  each  particle  and  their  behavior  is  based  on  an 
application  of  the  Betti-Rayleigh  theorem  from  classical  linear  elasticity  theory. 
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The  reciprocal  identity  or  the  Betti-Rayleigh  theorem  represents  a  relation  between  two 
different  static  elastic  states  on  the  same  body.  Consider  two  states,  both  in  a  state  of  static 
equilibrium,  (aj,  u/)  and  (a*1,  u/1)  with  surface  tractions 
t/  =  a;’  nj  and  t/1  =  a”  nj,  respectively.  Then  the  Betti-Rayleigh  theorem  states  that 

J  p  f'  u,11  dv  +  J  t/  uj11  ds  =  f  p  fj11  u/  dv  +  J  tin  u/  ds 

s  (H.1) 

where  v  denotes  the  volume  of  the  body  and  S  its  surface. 


Now,  from  virtual  work 

Jv  aij  eij!  dv  =  Jv  P  fi  dv  +  Js  ti  Uj 11  ds 

where  CTy  is  a  stress  state  in  equilibrium  and  u/1  is  a  displacement  field  and 
eij!  =  j  (uij  +  Uj1!).  Note  that  type  I  information  is  unrelated  to  type  II  information. 


01.2) 


Now  let  the  stress  state  I  be  of  the  form 


°ij  =  ^ii  ^lj 


The  corresponding  strains  are 

e  .i.  =  J-  ai _ 1 

1J  2  il  1J  2p(3X  +  2\x) 
Substitution  of  (II.3)  into  (II.4) 


I  _  1  51  S,  ^ 

ij  “  2(i  11  lj  ~  2p(3A.  +  2 p) 


01.3) 


^kk 


01-4) 


$kl  ^lk  Sjj 


01-5) 


One  displacement  field  which  satisfies  Eq.  015)  is 

X 


1  —  1  x  _ 

1  '  2 4  11  X|  2|1(32.  +  2(l) 


5;;  X; 


01.6) 


Now  by  using  Betti-Rayleigh  reciprocity  theorem  without  any  body  force  terms  f f  or  f/1, 
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Jv  °ij  e"  dv  =  |s  t/1  ui  ds 

Using  the  stress  state  in  Eq.  (II.  3) 

Jv  5ii  S]j  e"  dv  =  Js  t/1  u /  ds 
substitution  of  the  expression  of  u/  in  Eq.  (II.6)  gives 


Jy  Si,  8jj  Gjj1  dv  =  J  tj11 


1  s 

—  0U  Xj  - 
2(i 


- 7 -  5jj  Xj 

2(i(3X,  +  2(i)  1J  J 


ds 


Taking  type  II  information  to  be  the  change  in  the  actual  first  strain' rate  mode 


J,  <41  dv  =  L  <i 


1  s 

2?  5il  Xi  - 


2|i(3X  +  2n)  5il  Xj 


ds 


where  the  (•)  denotes  the  rate  of  change. 


Using  the  stress  state 


o\  =  Sj,  52j 


a  similar  expression  holds  for  the  second  strain  rate  mode 


JvE22<1v  =  Jtj 


1 


2(1  5i2  Xi  2(1  (3X  +  2(i)  ~1J  ~J 


8;;  X: 


ds 


01.7) 


01.8) 


01.9) 


(11.10) 


01.11) 


(11.12) 


The  stress  state 

Gjj  =  5U  52j 

gives 


^  2(i  5il52j  2(i  (3X  +  2(i)  §kl  §2k  5ij 


01.13) 


(11.14) 


The  second  term  in  the  right  hand  side  expression  is  identically  equal  to  zero,  therefore 
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(11.15) 


4  -  ^  5ii  h 


Using  a  displacement  associated  with  £jj  of  the  form 

u‘  =  ^  (Sii  52j  xj  +  Sji  52i  xj)  ai  16) 

the  third  strain  rate  mode  follows  using  the  same  argument 

Jv  ei2  dv  =  Js  i  (  8ii  xj  +  Sy  Xj  )  ds 

The  three  expressions  for  the  strain  rate  modes  all  may  be  written  in  a  single  expression 


fv  k  dv  =  Is [  (ti  xi  +  Xi)  " 


2ji(3X  +  4  n) 


tk  xk  djj]  ds 


01.18) 


For  the  particular  case  of  concentrated  (contact)  forces  applied  on  the  boundary 
ti  =  Fi  5(x  -  p)  8(y  -  ti)  8(z  -  y) 

Then  the  surface  integral  is  reduced  to  a  summation  over  the  contact  points  and 


(1119) 


kdv  =  X 

Contacts 


1 

4p 


(Fi  Xj  +  Fj  Xj)  + 


2p(3X  +  2  fi) 


Fk  xk  5ij 


(11.20) 


Equations  (11.20)  represent  three  constant  strain  rate  modes  in  each  grain.  These  constant 
strain  modes  will  be  used  to  update  the  comers  positions  as  well  as  the  internal  stresses  using 
Hooke’s  law. 


n.2.2  Evaluation  of  Force-Displacement  Relations  at  the  Contacts 

Local  deformations  are  assumed  to  occur  at  contacts  between  comers  and  edges,  see 
Fig.  II.l.  The  force-displacement  laws  are  evaluated  in  incremental  form,  that  is  to  say, 
changes  in  displacements  for  a  given  contact  are  evaluated  given  changes  in  horizontal,  and 
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vertical  displacements  and  rotations,  as  well  as  the  deformations  within  the  two  elements  in 
contact.  The  new  shear  and  normal  forces  are  then  calculated  by  summing  the  old  forces  and 
the  new  force  increments.  These  contact  forces  are  then  resolved  into  equivalent  horizontal 
and  vertical  forces  and  moments,  and  added  to  the  other  forces  and  moments  acting  on  the 
element.  In  the  original  PROBS  program  only  simplistic  linear  normal  contact  forces  were 
employed,  a  serious  limitation.  In  the  present  work  all  contact  forces  are  based  on  Mindlin’s 
theory  of  contact  [12]. 


eap  :  strain  rate 

Xq  :  global  centroid  coordinate  of  element  i 

Xi  :  global  centroid  coordinate  of  element  j 

X£  :  global  coordinate  of  contact  point 


The  incremental  shear  and  normal  displacements  are  therefore 
AXa  =  Jap  Xp  At 

where 


(11.22) 
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[J«P]  = 


cosy  smy 
-siny  cosy 


01.23) 


and  y  is  the  angle  of  the  edge  with  the  global  X  coordinate  axis. 


The  contact  shear  T  and  normal  force  N  at  the  contact  are  based  on  Mindlin’s  theory 
[12]  and  are  evaluated  assuming  rounded  contact  surfaces 


N  = 


4RiRj 


1  -v2 


9  nr  (Ri  +  Rj)  j 


1/2 


54/2 


(II.24a) 


T  = 


-  f  1 - >2 

l  KfN^3  j 


KfN 


f  N  (  |  T  |  <  f  N  ) 


(II.24b) 


or 


|  T  |  =  fN 


(  !T|  >  fN) 


(n.24c) 


In  these  formulas  K,  R,  and  Rj  are  given  by  the  expressions 


3(2  -  v) 


IV  - 

16(i  ( 

6n2E  RjRj 

(11.25 a) 

■  8(1  -  V2)(R;  +  Rj)  J 

Ri  = 

„  K 

j 1/2  =  Average  Radius  of  Element  i 

(n.25b) 

Rj  = 

'A 

.  7t 

j1/2  =  Average  Radius  of  Element  j 

(II.25c) 

and  Aj,  Aj  are  the  areas  of  elements  i  and  j,  respectively. 

coefficient  between  the  contacting  bodies. 

Also,  f  denotes  the  friction 
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n.2.3  Numerical  Integration  of  the  Equations  of  Motion 

The  force  and  moments  generated  on  each  element  are  used  to  compute  horizontal, 
vertical  and  rotational  accelerations.  The  accelerations  are  then  integrated  to  give  velocities, 
and  integrated  again  to  give  displacements. 


Each  grain,  with  mass  n^  and  moment  of  inertia  I;,  is  assigned  three  independent 
degrees  of  freedom  at  the  element  centroid,  namely  Xlt  X2  and  0,  the  displacements  in  the  1 
and  2  directions  and  rotation,  respectively.  The  grain  must  obey  the  following  equations  of 
motion 


mi 


dzX0 

dt2 


+  q 


dXo 

dt 


+  Fia  =  0  a=l,2 


(11.26) 


d20 

dt2 


+  Di^ 

1  dt 


+  Mi  =  0 


where  Fia  is  the  sum  of  forces  in  the  1  or  2  direction,  and  M,  is  the  sum  of  all  moments  of 
force  (about  the  centroid)  for  the  particle.  The  constants  Q  and  Dj  denote  damping 
coefficients. 


Consider  the  nonlinear  second  order  differential  equation  a  special  case  of  which  is 
any  of  Eq.  (11.26), 

X  +  a  X  =  —  F(X,  X) 

m  (11.27) 

In  Eq.  (11.27)  X,  X,  X  are  displacement,  velocity  and  acceleration,  respectively,  and  m  is  the 
mass  of  the  element.  The  quantity  F(X,  X)  is  the  force  generated  from  contacts  and  other 
applied  (body)  forces. 


In  order  to  integrate  the  equations  the  derivatives  at  time  t  may  be  written 
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A  +  At 


X  = 


X  = 


-  X1 


At 


Xl  +  At  +  X1 


(II.28a) 


(II.28b) 


Substitution  of  Eqs.  (II.28a)  into  Eq.  (II.28b)  gives  the  following 

Xt  +  At- 
At 


Cl  f  Xt  +  At  +  xc  1 

1 

-  +  a  - - - 

=  —  F(X,X) 

l  2  < 

m 

velocity  at  t  +  At: 

- 

f  1  -  a  At/2  ]  •  t 

r  At  1  F(X,  X) 

l  1  +  a  At/2  J 

.  1  +  a  At/2  J  m 

(11.29) 


(11.30) 

The  new  values  for  velocities  are  used  to  update  the  displacements  and  rotations  of  the 
elements  by  a  further  numerical  integration 

Xt  +  At  =  X1  +  XAt 

(0.31) 


This  integration  scheme  is  unconditionally  stable  if  F  is  independent  of  X  or  X.  If  the 
force  is  a  linear  function  of  X,  i.e.,  linear  elastic  behavior  is  present  (F  =  KX),  then  the  pro¬ 
cedure  is  conditionally  stable  provided  At  <  2Vm/K.  But  the  presence  of  the  nonlinear  force 
term  F(X,X),  which  depends  strongly  on  the  number  of  contacts  for  each  grain,  may  produce 
numerical  instability  in  the  solution.  Under  numerically  unstable  conditions  the  contact  forces 
will  be  very  large,  and  fluctuate  wildly  from  iteration  to  iteration.  (The  results  generated  by 
the  program  are  meaningless  under  such  a  situation;  the  problem  must  be  re-run  with  a 
smaller  time  step.) 


n.2.4  Rotation  Effects  on  Stress  Calculations  in  Each  Grain 


If  a  grain,  together  with  the  forces  acting  on  it,  is  rotated,  the  internal  stresses,  referred 
to  local  axes,  are  unaffected.  However,  the  stresses,  as  expressed  in  global  coordinates  will 
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change.  It  is  common  in  Lagrangian  finite-difference  codes  to  apply  stress  correction  terms  to 
the  stresses  at  each  time  step  in  order  to  allow  for  this  apparent  stress  change  due  to  rotation. 
The  correction  terms  are  derived  from  the  tensor  transformation  equations  used  to  determine 
stresses  in  a  new  coordinate  system  when  they  are  known  in  an  old  system. 


Denote  the  old  coordinates  system  by  X;  (i=l  ,2,3)  and  the  new  coordinates  system  by 
Xj  (i=l,2,3).  The  tensorial  relations  that  relate  the  stresses  cfy  in  the  new  coordinate  and  the 
stresses  Oy  in  the  old  coordinate  system  are  as  follows. 

ax'; 

°'J  ~  °aP  d*a  dXp  01.32) 


Let 


cos9  sin9 
-sin9  cos0  . 

where  0  is  the  angle  of  rotation  between  the  two  coordinate  systems, 
say  A0  occurring  over  a  single  time  step  At,  [Jy]  can  be  approximated 


3X: 


01.33) 

If  0  is  a  small  change, 
by 


[Jij] 


1  A0 
.  -A0  1 


01.34) 


Noting  that  A0  =  0  At,  then 


Jij  =  8jj  +  Ry  At 

where 


01.35) 


[Rij]  = 


0  0 
— 0  0 


01.36) 


Using  this  expression  for  Jy  in  (11.35) 


19 


®ij  ^a|J  (^ia  ^-ia  At)(8jp  +  Rj|J  At) 

—  ®ij  (^aj  Rja  ■*"  ^ip  ^jj})  ^ 


(11.37) 

These  new  stresses  correspond  to  a  change  in  coordinate  axes.  However,  in  the  case  of  a  par¬ 
ticle  rotating  and  the  coordinate  axes  remain  fixed,  and  the  correction  term  is  opposite 


AOjj  -  -  (<Jaj  Rja  +  CJjp  Rjp)  At 


(11.38) 


In  MASOM  the  internal  stresses  are  updated  in  each  element  using  the  strain  rate  and 
stress  correction  as  follows 


=  Ojj  +  (  X  Efck  8y  +  2(j.  £y  -  RiK  -  CTjlK  RiK  )  At 


(11.39) 


II.2.5  Fracture  of  the  Elements  in  MASOM 


The  fracture  of  grains  under  multiaxial  loading  conditions  can  be  an  important  aspect 
of  the  micro-mechanics  of  frictional  materials.  The  grain  fracture  criterion  incorporated  in 
MASOM  is  based  on  the  assumption  that  cracking  occurs  when  point  loads  are  applied  on 
opposite  sides  of  a  grain.  In  its  ideal  form,  this  type  of  approach  constitutes  the  "Brazilian 
test"  and  is  performed  with  a  disc  loaded  between  flat  plates.  More  generally,  many  point 
load  tests  have  been  performed  in  which  irregular  shapes  are  loaded  between  two  identors  that 
may  have  small  radii  of  curvature.  These  tests  are  characterized  by  a  failure  load  that  is 
related  to  the  distance  between  the  points  of  application  of  the  loads,  and  a  strength  that  is 
supposed  to  be  constant  for  a  given  material. 


In  the  MASOM  program  the  relationship  for  fracture  is  taken  to  be  the  simple  form 
S  =  K  F/d2,  where  F  is  the  average  of  co-linear  compressive  forces  applied  on  opposite  sides 
of  the  grain,  d  is  the  distance  between  applied  loads,  k  is  a  factor  that  depends  on  the  shape 
of  the  grain,  and  S  is  the  material  strength,  taken  to  be  constant. 


During  each  cycle  of  the  solution,  each  grain  is  scanned  for  the  sum  of  the  largest  co- 
linear  forces,  say  Ft  and  F2  applied  to  the  grain  by  comers  from  other  grains.  The  average  of 
these  forces  xh  (F;  +  F2)/2  is  then  taken,  divided  by  their  separation  distance  (squared),  and 
compared  to  the  specified  strength  S.  If  Vi  (Fj  +  F2)/d2  >  S  crack  is  formed;  if 

V2  (Fj  +  F2)/d2  <  S  no  crack  is  formed.  Once  the  crack  has  been  formed;  the  element 

is  separated  into  two  pieces  along  the  line  of  the  crack.  The  coordinates  of  the  new  elements 
are  reset  and  the  new  grains  are  created.  Internal  stress  and  strains  are  set  equal  zero  in  the 
new  grains. 

II.3  Analysis  of  Material  Behavior  Using  MASOM 
II  .3.1  Introduction 

The  computer  program  MASOM  models  the  movement  and  interaction  of  the  many 
particles  which  make  up  the  material,  during  simulations  of  material  response.  The  sample  of 
the  material  is  generated  from  a  microphotograph  of  a  real  granular  material;  each  particle  is 
modeled  as  a  polygon  of  a  shape  and  size  resembling  a  particle  in  the  photograph.  A  digitizer 
is  used  to  generate  the  data  file  to  be  used  by  the  computer  program. 

During  such  numerical  experiments,  a  variety  of  boundary  conditions  and  loading  as 
well  as  unloading  sequences  may  be  specified,  and  measurements  made  of  contact  force  orien¬ 
tation  and  magnitude,  as  well  as  average  stress  and  strain,  or  other  variables  of  interest. 

In  this  section  the  physical  tests  which  simulated  using  MASOM  are  the  uniaxial  strain 
and  triaxial  compression  and  simple  shear  tests. 
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H.3.2  Uniaxiai  Test  Without  Fracture 

The  uniaxial  test  simulation  was  based  on  data  for  Cambria  beach  sand.  The  computer 
model  in  Fig.  II.2  consists  of  227  particles  with  maximum  dimension  varying  between  0.5  and 
2.0  mm.  The  elastic  properties  for  the  grains  are  [37], 
E  =  4.5  x  107  Kg/mm  •  s2.  Young’s  modulus,  and  G  =  1.8  x  107  Kg/mm  •  s2,  the 
Shear  modulus.  Three  additional  elements  are  fixed  to  simulate  the  fixed  uniaxial  strain  test 
chamber  and  two  carry  the  loads  which  are  applied  gradually  over  time.  Figure  n.3  shows 
the  sample  with  a  gravity  field  acting  on  it  but  without  any  load  applied  by  the  apparatus. 
The  intensity  of  the  shading  gives  an  indication  on  the  magnitude  of  the  different  stresses 
inside  each  grain.  (Vertical  shading  correspond  to  the  vertical  stress  oy,  horizontal  shading 
correspond  to  cx,  and  the  oblique  shading  correspond  to  the  shear  stress  crxy.)  Note  the  build 
up  of  stress  with  depth,  i.e.,  the  lithostatic  effect.  Figure  II.3  shows  concentrations  of  stiff 
chains  of  particles  that  carry  the  applied  load  (in  this  case  is  the  weight  of  the  particles), 
which  is  very  non  uniform.  Figure  II.4a  and  II.4b  show  small  dots  on  some  of  the  comers, 
the  dots  indicating  that  the  corresponding  comers  are  sliding  along  their  contacting  edges. 
Evidently  lithostatic  behavior  is  highly  inhomogeneous. 

The  uniaxial  strain  test  was  simulated  by  integrating  the  system  for  10,000  cycles. 
Contact  force  magnitudes  and  directions  as  well  as  displacements  and  velocities  were  col¬ 
lected  every  100  cycles  to  construct  the  stress  strain  relation  for  the  test.  Vertical  stress 
versus  vertical  strain  curve  is  shown  in  Fig.  II.5.  The  maximum  strain  reached  about  11% 
and  the  residual  strain  obtained  after  unloading  was  around  3.5%.  Figure  11.6  shows  the  stress 

ax  +  2  ay 

path  which  is  the  stress  difference  or  deviatoric  stress  ox  -  cy  versus  "pressure"  - - - . 

Since  the  program  is  2D  there  is  a  tacit  assumption  that  ay  =  az,  and  this  assumption  is  used 
whenever  pressure  effects  are  quantified. 
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Figure  II.  7  shows  the  variation  of  the  number  of  contacts  as  a  function  of  the  vertical 
stress.  The  relationship  between  stress  and  the  number  of  the  contacts  is,  with  the  exception 
of  early  redistributions,  remarkably  close  to  that  between  stress  and  strain. 

IL3.3  Uniaxial  Test  With  Fracture 

The  uniaxial  test  simulation  with  fracture  allowed  was  re-analysed  using  the  same 
Cambria  beach  sand  material  properties  as  in  H.3.2.  In  addition  to  these  elastic  properties,  the 
material  strength  constant  was  set  S  =  108  Kg/mm  s2.  Recognizing  that  the  cracking  limit  is 
inversely  proportional  to  the  square  of  the  distance  between  the  application  of  the  contact 
forces,  this  criterion  results  in  fracture  of  the  smaller  particles  first.  For  this  uniaxial  test  47 
particles  have  fractured  out  of  227  particles  for  the  model.  The  maximum  cracking  line  was 
1.15  mm  (this  model  the  maximum  dimension  for  the  particles  size  varies  between  0.5  mm 
and  2.0  mm). 

The  plots  shown  in  Figs.  II. 9  through  11.14  display  results  for  the  uniaxial  test.  Frac¬ 
ture  occurs  around  the  cycle  1400  when  the  material  fracture  strength  is  reached.  Initiation  of 
cracking  occurs  first  on  the  smaller  particles,  which  because  of  their  small  dimension  are 
unable  to  resist  intense  loads  as  effectively  larger  grains  with  equal  strength  properties.  As 
cracking  occurs  some  re-distribution  of  forces  develop  but  the  general  patterns  of  contact 
force  and  force  direction  are  disturbed  only  slightly.  The  unloading  phase  begins  after  cycle 
1600. 


The  plots  in  Figs.  11.15  through  11.17  displays  transmission  of  the  contact  forces  inside 
each  element.  The  thickness  of  these  lines  are  proportional  to  the  magnitudes  of  the  contact 
forces.  A  remarkable  feature  here  is  the  preponderance  of  vertically  oriented  contact  forces, 
i.e.,  the  low  material  confinement  forces. 
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Contact  force  magnitudes  and  directions  as  well  as  displacements  and  velocities  were 
collected  for  all  particles  of  the  sample  and  at  the  boundaries.  These  forces  and  displace¬ 
ments  are  used  to  construct  an  overall  stress  strain  relation  for  the  material.  Vertical  stress 
versus  vertical  strain  is  shown  in  Fig.  11.18.  The  curve  shows  a  sudden  drop  of  the  stress 
level  near  the  peak,  i.e.,  when  the  particles  are  allowed  to  fracture.  After  rearrangement  of 
the  particles  due  to  fracture,  the  sample  will  regain  some  strength  before  the  unloading  phase. 
The  maximum  strain  reached  was  about  18%  and  the  residual  strain  obtained  after  unloading 
was  about  14%.  Figure  0.19  shows  the  stress  path  taken,  i.e.,  the -stress  difference  or  devia¬ 
te  +  2ay 

toric  stress  (ox  -  cy)  versus  pressure  ( - - — -). 

IL3.4  Qualitative  Comparison  of  MASOM  Results  With  Laboratory  Tests 

Figure  11.20  shows  the  stress  strain  relation  obtained  in  the  UCLA  soils  laboratory  for 
the  uniaxial  strain  test  on  Cambria  Sand.  The  test  was  performed  using  a  tri axial  cell  and  by 
applying  gradually  increasing  lateral  stresses  in  order  to  maintain  zero  lateral  strain.  In  spite 
of  differences  between  the  void  ratio  in  the  laboratory  sample  and  the  computer  model  and  in 
applied  stress  levels  (150  bar  maximum  stress  in  the  computer  simulation  compared  to  25  bar 
in  the  laboratory  test),  the  numerical  simulation  appears  to  have  captured  the  general  behavior 
of  the  material  for  the  loading  phase  as  well  as  the  unloading  phase,  see  Fig.  11.18  and  Fig. 
11.20. 

II.3.5  Triaxial  Compression  Test 

The  triaxial  compression  test  using  a  sample  of  Cambria  sand  was  analysed  using 
MASOM  but  without  fracture.  The  computer  model  in  Fig.  11.21  illustrates  the  boundary  con¬ 
ditions  used  for  the  numerical  test.  The  horizontal  plates  carry  the  vertical  loads,  both  top  and 
bottom  are  allowed  to  move.  Confining  loads  are  carried  by  the  vertical  elements. 
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In  the  first  phase  of  loading,  that  of  isotropic  compression,  vertical  and  horizontal 
loads  are  increased  gradually  and  at  the  same  rate  up  to  50  bar.  In  the  second  phase,  termed 
triaxial  compression,  the  vertical  force  is  increased  up  to  150  bar  while  maintaining  the  same 
confining  pressure  of  50  bar  in  the  horizontal  direction.  The  increasing  vertical  load  produces 
a  downward  motion  of  the  top  plate,  but  variable  behavior  of  the  vertical  confinement. 

The  plots  shown  in  Fig.  11.22  through  11.32  displays  results  for  both  phases  of  the 
triaxial  compression  test  with  the  stress  level  reached  in  both  the  vertical  and  horizontal  direc¬ 
tions.  Figures  n.33  through  11.39  show  the  contact  forces  transmission  displayed  inside  each 
particle.  The  thicknesses  of  these  lines  are  proportional  to  the  magnitudes  of  the  contact 
forces. 


The  element  data  and  contact  forces  information  were  collected  during  the  numerical 
analysis.  Contact  forces  and  boundary  displacements  were  used  to  construct  the  curve  in  Fig. 
11.40  which  represents  the  deviatoric  stress  (a}  -  c2)  versus  the  vertical  strain  Figure  11.41 
shows  the  stress  path  versus  V2  •  a2.  Clearly  the  first  section  of  the  curve  represents  the 
isotropic  compression  phase  and  the  second  section  represents  the  vertical  loading  phase. 

II.3.6  Observed  Behavior  from  the  Numerical  Test 

The  computer  simulation  using  MASOM  is  shown  to  realistically  capture  macroscopic 
behavior  observed  in  physical  tests.  The  effects  of  several  important  processes  are  revealed 
by  the  numerical  simulation,  namely 

1.  The  effect  of  particle  movement  (displacement  and  rotation). 

2.  The  particles  elastic  behavior  and  subsequent  fracture  at  high  stress  levels. 

3.  Inhomogenous  behavior  of  the  system  of  particles  during  loading. 

4.  Build  up  of  loading  chains  and  islands  of  quasistatic  subassemblies. 
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5.  Relationship  between  intergranular  contact  numbers  and  strain  resembles  the  stress 

strain  relationship. 

n.3.7  Evolution  of  Particulate  Mass  During  Loading 

The  behavior  of  the  frictional  material  is  very  complex,  involving  a  number  of 
different  physical  processes  as  may  be  seen  from  the  computer  generated  field  plots  of  the 
particulate  assembly  at  various  phases  of  loading.  To  gain  a  better  understanding  of  the  evo¬ 
lution  of  particle  interaction,  a  video  tape  was  generated  of  the  response  of  the  system 
throughout  the  duration  of  the  load-unload  cycle.  The  evolution  of  physics  observed  during 
the  numerical  experiment  shows  that  the  applied  loads  are  carried  by  chains  of  grains  in  the 
direction  of  maximum  compressive  strain  rate.  The  distortion  of  these  chains  resembles  the 
buckling  of  a  curved  beam  or  arch,  and  neighboring  particles  exhibit  large  motion  and  mainly 
rotation  without  the  appearance  of  large  stresses. 

The  structural  behavior  of  the  granular  material  considered  appears  to  consist  of  two 
structures  behaving  differently.  The  first  is  made  up  of  stiff  chains  and  behaves  like  a  com¬ 
plex  elastic  truss.  The  deformation  patterns  of  this  truss  imposes  local  deformation  on  the 
particles  of  a  second  type  of  structure  confined  by  the  chains  of  the  first  structure  whose 
behavior  is  similar  to  a  perfectly  plastic  material  due  to  the  relatively  small  compressive  stress 
acting  on  it.  As  the  applied  load  increases  the  elements  of  the  first  structure  become  progres¬ 
sively  more  kinked,  and  therefore  less  able  to  carry  load.  The  load  is  then  transferred  to 
neighboring  elements  of  the  secondary  structure  as  increasing  deformation  brings  more  grains 
into  contact. 

II.4  Concluding  Remarks 

The  micromechanical  analysis  conducted  using  the  MASOM  computer  program  reveals 
an  intricate,  and  continually  changing  microphysics  at  the  intra-  and  interparticle  level.  The 
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program  permits  a  careful  examination  of  granular  displacement  and  rotation,  deformation  and 
possible  fracture  at  high  stress  levels.  Given  the  results  of  the  numerical  simulation,  which 
includes  contact  force  distribution  and  displacements  along  the  boundaries,  it  is  possible  to 
construct  the  stress  strain  relations  for  overall  behavior  of  finite  samples  of  the  frictional 
material.  In  particular,  two  standard  laboratory  tests,  the  uniaxial  strain  test  and  the  triaxial 
compression  test  were  analysed  numerically  for  loadings  well  into  the  inelastic  range.  The 
numerical  analysis  has  captured  the  general  behavior  of  the  material  during  both  loading  and 
unloading  phases. 

However,  it  is  clear  that  the  present  analysis  could  be  improved  in  future  work. 
Several  areas  uf  research  should  be  mentioned  where  real  progress  could  be  made  to  better 
understand  the  micro-mechanical  behavior  of  granular  materials.  First  and  simplest,  the 
current  sample  used  in  both  the  uniaxial  strain  test  and  the  triaxial  compression  test,  consist¬ 
ing  of  227  particles,  is  too  small  to  permit  the  investigation  of  stability  issues.  In  the  current 
research,  the  size  of  the  computer  model  was  restricted  due  to  limited  computational  resources 
and  the  numerically  intensive  calculations  involved  for  the  type  of  simulation  using  MASOM. 
It  is  hoped  that  more  computational  support  will  be  available  in  the  future  to  permit  more 
finely  resolved  models  with  more  grains,  so  that  stability  issues  can  be  explored  numerically. 

Another  important  computer  resource  constraint  pertains  to  the  modeling  of  individual 
grains.  In  the  present  model  the  polygonal  shapes  were  restricted  to  a  maximum  of  10  edges 
per  particle.  A  smoother  representation  of  the  boundary  of  the  particle  may  significantly 
affect  the  response  of  the  system  especially  at  high  stress  levels  and  for  cases  involving  angu¬ 
lar  particles  with  relatively  low  fracture  strengths. 

Other  more  advanced  physics  can  also  be  implemented  in  MASOM  which  will  be 
helpful  to  better  understanding  material  behavior.  For  example: 
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1. 


Specification  of  pressure  profile  for  the  boundary  conditions.  At  this  time  only  point 
loads  can  be  applied  to  the  boundaries.  This  limits  the  shape  and  form  of  boundary 
support,  a  serious  restriction,  especially  with  regard  to  types  of  failures  where  fracture 
planes  or  zones  extend  to  the  boundaries  of  the  sample,  i.e.,  interact  with  the  boundary 
conditions. 

2.  Motion  of  fluid  flow  through  the  joints  and  voids.  Two  phase  flow  can  be  modelled  by 
combining  pressure-flow  calculation  with  force-displacements  calculation.  This 
research  has  already  begun. 

3.  Extension  of  MASOM  from  2D  to  3D.  Of  course,  the  real  world  is  three  dimensional 
and  it  would  be  very  helpful  if  a  discrete  element  program  could  be  developed  for 
investigating  real  particles,  and  more  important,  3D  particulate  assemblies  where  inter¬ 
stitial  behavior  (filing  from  crushed  particles,  flow  to  fluids,  etc.)  is  probably 
significantly  different  from  a  2D  model.  The  primary  difficulty  is  simply  managing 
(i.e.,  minimizing)  the  computational  burden  as  polygons  become  polyheldrons  and  the 
contact  problem  becomes  extraordinarily  complex. 
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II. 2  Uniaxial  strain  test  Computer  model. 
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Uniaxial  strain  test.  Sliding  of  particles,  cycle  800 
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II.4b 


Uniaxial  strain  test.  Sliding  of  particles,  cycle  1200 
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6  Uniaxial  test:  Stress  path,  deviatoric  stress  versus  pressure. 
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Uniaxial  test.  Vertical  stress  versus  total  number  of  contacts. 


Uniaxial  test  with  fracture.  Computer  model. 
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II- 9  Uniaxial  test  simulation,  0  bar  <  a£  <  1  bar 
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II.  10  Uniaxial  test  simulation,  2  bar  <  <  10  bar 
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11.12  Uniaxial  test  simulation,  90  bar  <  Oj  <  150  bar 
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11.13  Uniaxial  test  simulation,  90  bar  <  (jj  <  120  bar 
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11.15  Uniaxial  test  simulation.  Contact  force  display 
cycle  1920-2010 
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11.16  Uniaxial  test  simulation.  Contact  force  display 
cycle  2040-2130 
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11.17  Uniaxial  test  simulation.  Contact  force  display 
cycle  2160-2250 
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Uniaxial  test  with  fracture.  Deviatoric  stress  versus  pressure. 
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n.21  Triaxial  test  simulation.  Computer  model. 
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11.22  Triaxial  test  simulation,  0  bar  <  (*i  <  5  bar,  0  bar  <  ct2  <  5  bar 
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n.23  Triaxial  test  simulation,  5  bar  <  <  10  bar,  5  bar  <  <  10 
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11.24  Triaxial  test  simulation,  10  bar  <  <  20  bar,  10  bar  <  o2  <  20 
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11.25  Triaxial  test  simulation,  20  bar  <  <  30  bar,  20  bar  <  a2  <  30 

bar 
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11.26  Triaxial  test  simulation,  30  bar  <  Oj  <  50  bar,  30  bar  <  02  <  50 
bar 
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Triaxial  test  simulation,  70  bar  <  cii  <  90  bar,  c2  =  50  bar 
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11.29  Triaxial  test  simulation,  90  bar  <  Gj  <  110  bar,  cr2  =  50  bar 
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11.32  Triaxial  test  simulation,  120  bar  <  (Ji  <  130  bar,  ct2  =  50  bar 
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11.33  Triaxial  test  simulation.  Contact  force  display 
cycle  1020-1080 
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Triaxial  test  simulation.  Contact  force  display 
cycle  1100-1160 
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Triaxial  test  simulation.  Contact  force  display 
cycle  1180-1240 
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Triaxial  test  simulation.  Contact  force  display 
cycle  1340-1400 
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0.38  Triaxial  test  simulation.  Contact  force  display 
cycle  1420-1480 
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11.39  Triaxial  test  simulation.  Contact  force  display 
|  cycle  1500-1560 
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Triaxial  test. 

Deviatoric  stress  versus  vertical  strain. 
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11.41  Triaxial  test.  Vertical  stress  versus  confining  pressure. 
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HI.  DEVELOPMENT  OF  THE  EXPERIMENTAL  EQUIPMENT 


HI.1  Introduction 

The  experimental  equipment  for  triaxial  testing  at  high  pressures  has  been  assembled 
and  is  suitable  to  performance  of  triaxial  tests  as  originally  envisioned.  The  equipment  con¬ 
sists  of  a  large  frame  with  a  hydraulic  cylinder  for  vertical  loading,  and  pressure  supply  units 
for  lateral  as  well  as  vertical  loading.  An  integral  portion  of  the  testing  equipment  consists  of 
a  dedicated  computer  which  is  used  for  datalogging  and  for  control  of  the  experiments.  Any 
conceivable  stress-path  or  strain-path  possible  in  the  triaxial  apparatus  can  be  followed  under 
stress  control  or  strain  control.  Figure  1  shows  a  schematic  diagram  of  the  experimental 
setup. 


As  described  in  previous  progress  reports,  the  first  as  well  as  the  second  design  of  a 
loading  frame  both  turned  out  to  be  prohibitive  in  terms  of  cost.  Alternative  solutions  were 
sought.  Several  used  loading  frames  came  on  the  market  in  the  Spring  of  1987,  and  one  of 
these  fulfilled  the  requirements  for  this  project.  This  loading  frame  was  first  inspected  (by 
Lade)  in  June  at  General  Electric  in  Schenectady,  N.Y.,  and  subsequently  purchased  and  tran¬ 
sported  to  UCLA  were  it  arrived  in  August  of  1987. 

The  loading  machine  was  fabricated  by  Gilmore  (in  Cleveland,  Ohio)  and  purchased  as 
one  of  three  by  General  Electric  in  1972.  It  has  a  capacity  of  200  kips,  and  it  is  relatively 
stiff,  an  important  requirement  for  controlled  descent  along  the  strain  softening  branch  of  any 
force-deflection  curve.  This  loading  frame,  which  came  with  a  hydraulic  cylinder  for  vertical 
loading,  was  substantially  cheaper  than  any  frame  that  could  be  designed  and  produced  for 
this  project.  Further,  the  machine  does  not  appear  to  have  been  used  heavily,  and  all  parts 
have  been  checked  and  are  in  good  working  condition. 
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While  the  loading  frame  was  being  acquired,  the  equipment  to  produce  pressurized  oil 
and  cell  fluid  (water)  was  purchased.  The  pressure  required  in  the  hydraulic  cylinder  as  well 
as  the  confining  pressure  in  the  triaxial  cell  are  generated  by  step  motors  which  act  on  ball 
screws  through  gears  to  compress  oil  and  water  in  smaller  hydraulic  cylinders.  In  addition,  a 
pressure  booster  pump  to  produce  constant  high  pressure  in  water  (from  a  low  input  air  pres¬ 
sure)  may  be  used  for  conventional  triaxial  tests  performed  at  constant  pressure.  The  actually 
achievable  maximum  cell  pressure  is  15,000  psi,  as  desired. 

Both  the  vertical  load  and  the  triaxial  cell  pressure  are  controlled  by  the  step  motors  to 
achieve  any  desirable  stress  path  or  strain  path.  This  i«  accomplished  using  an  integrated  data 
acquisition  system  run  by  the  same  computer  that  controls  the  step  motors. 

The  test  data  is  collected  from  pressure  transducers,  LVDTs  and  load  cells.  In  addi¬ 
tion,  quantities  such  as  pressures  and  deformations  will  be  monitored  visually  on  pressure 
gages,  and  dial  gages.  All  high  pressure  lines  are  provided  with  relief  valves  for  safety. 

A  considerable  amount  of  work  including  purchasing  of  computer  hardware  and 
software,  as  well  as  programming  for  test  control,  data  logging  and  processing  has  been  per¬ 
formed  by  two  graduate  students.  The  high  levels  of  quality  and  sophistication  in  this  high 
pressure  testing  equipment  is  well  supported  by  the  acquisition  of  the  professionally  built,  but 
inexpensive  loading  machine. 

The  triaxial  cell  in  which  the  specimens  are  tested  has  been  constructed.  This  cell  is 
designed  for  a  maximum  cell  pressure  of  15,000  psi  (1000  bars).  The  maximum  vertical  load 
supplied  by  the  loading  machine  (200,000  lbs)  is  more  than  sufficient  to  produce  failure  in  a 
3.0  inch  diameter  triaxial  specimen  under  undrained  conditions  at  the  maximum  cell  pressure. 
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Each  of  the  major  components  of  the  experimental  equipment  is  described  in  further 
details  in  the  following  chapter.  This  includes  the  loading  frame,  the  high  pressure  triaxial 
apparatus,  the  pressure  sources,  and  the  test  control  panel.  The  computerized  datalogging  and 
test  control  system  is  explained  in  a  separate  chapter  below.  Finally  the  results  of  a  drained 
test  and  an  undrained  test  produced  by  the  newly  developed  equipment  are  presented  to 
demonstrate  its  capabilities. 

HI.2  Details  of  Experimental  Setup 

III.2.1  Description  of  the  Loading  Frame 

The  loading  frame  is  a  "Gilmore  model  426-200".  It  is  designed  to  withstand  200 
Kips  static  loads  and  100  Kips  dynamic  loads.  The  frame  consists  of  a  base,  four  guide 
columns,  and  a  crosshead. 

The  base  is  a  51"W  x  35"D  x  42"H  welded  steel  fabrication.  It  serves  as  an  anchor 
for  the  guide  columns  and  a  housing  for  the  actuator.  Four  guide  columns  are  anchored  in  the 
base.  Each  column  is  five  inches  in  diameter  and  is  of  sufficient  length  to  allow  the  test  zone 
to  be  adjusted  from  zero  to  77  inches.  Maximum  overall  height  of  the  columns,  including 
base,  is  161  inches.  The  crosshead  is  guided  on  the  four  columns.  It  is  a  52"W  x  33"D  x 
15"H  welded  steel  fabrication.  The  crosshead  is  clamped  to  the  columns  by  14  clamp  studs. 
It  is  lifted  into  position  by  two  hydraulic  cylinders  mounted  centrally  on  each  side  of  the 
machine. 

The  actuator  used  in  the  loading  frame  is  a  "Gilmore  model  L433A-180"  hydraulic 
actuator.  It  is  of  the  double  acting  type  with  a  double  ended  piston  rod.  The  rating  of  the 
actuator  is  200,000  lb  static  or  180,000  lb  dynamic  at  an  operating  pressure  of  3000  psi.  The 
total  stroke  of  the  piston  is  6  inches.  The  pressure  in  the  actuator  is  measured  by  a  "HBM" 
pressure  transducer  located  outside  the  piston  and  by  a  "Heise"  dial  gage  on  the  control  panel. 
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Fixed  on  the  crosshead  is  a  20  Ton  "HBM  model  C3H2"  load  cell  that  is  used  to 
measure  vertical  loads.  At  higher  loads,  a  50  Ton  "HBM  model  C3H2"  and  a  100  Ton  "Inter¬ 
face  model  1240-AF"  load  cells  are  used. 

Two  500  lb  "DUFFLYNX"  electrical  wire  rope  hoists  were  installed  to  lift  the  high 
pressure  triaxial  apparatus  up  on  the  pedestal  on  the  loading  machine.  The  hoists  are  mounted 
on  the  crosshead  and  slides  on  two  6  feet  long  I  beams. 

IIL3  Description  of  the  High  Pressure  Triaxial  Apparatus  (HPT- A) 

The  HPTA  when  assembled  stands  23  in.  high,  is  12  in.  in  diameter,  and  weights 
approximately  450  lb.  It  has  four  major  components:  the  upper  cap  plate,  the  loading  piston, 
the  base  plate,  and  the  cylinder  wall.  All  four  pieces  are  fabricated  from  steel  AISI  4340. 
They  were  heat  treated  to  a  yield  strength  of  180  Ksi  with  a  Rockwell  hardness  of  40/43. 

HI.3.1  Upper  Cap  Plate  and  Loading  Piston 

The  upper  cap  plate  houses  a  linear  bearing  "Boston  LBB-3000  extra-precision".  The 
linear  bearing  guides  a  17  in.  long  loading  piston  with  3  in.  diameter. 

The  piston  is  dense  chromed  with  a  surface  finish  of  4  RMS  to  provide  a  sufficiently 
smooth  surface  against  which  to  seal.  A  reciprocal  seal  made  by  "Ball  Seal  Engineering"  is 
used  to  seal  15,000  psi  of  fluid  pressure  inside  the  triaxial  cell.  It  is  made  from  graphite 
fiber-reinforced  PTFE  and  is  activated  by  a  medium-load  spring.  A  backup  ring  made  from 
high  strength  polymer  is  placed  behind  the  seal  to  reduce  extrusion  during  high  pressure 
operation.  The  seal  is  held  in  place  by  an  aluminum  ring  attached  to  the  upper  plate  by  6 
screws.  It  can  be  removed  for  replacement. 
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The  upper  cap  plate  has  one  groove  that  houses  an  O-ring  which  seals  the  fluid  inside 
the  cylinder  wall. 

The  top  cap  incorporates  a  single  drainage  line  used  to  vent  any  air  trapped  inside  the 
cell  prior  to  testing.  The  drainage  line  is  closed  by  a  high  pressure  valve  located  at  its  port 
outside  the  cap. 

ni.3.2  The  Cylinder  Wall 

The  cylinder  wall  has  an  inner  diameter  of  5  in.  and  a  thickness  of  1  in.  The  triaxial 
cell  is  designed  for  a  pressure  of  15,000  psi  with  a  factor  of  safety  of  2.5. 

m.3.3  The  Base  Plate 

The  lower  base  plate  incorporates  a  single  chamber  fluid  line  and  two  specimen 
drainage  lines.  The  fluid  line  is  used  to  pressurize  the  cell,  and  the  fluid  pressure  is  measured 
by  an  "HBM"  pressure  transducer  installed  at  the  outside  port  of  the  line.  The  two  drainage 
lines  are  connected  to  the  soil  specimen  top  and  bottom  and  plates.  When  opened,  they  con¬ 
duct  fluid  from  the  specimen  to  the  volume  change  device.  A  pressure  transducer  located  on 
the  drainage  line  outside  the  cell  is  used  to  measure  the  pore  pressure  in  undrained  tests. 

The  lower  base  plate  has  three  grooves  machined  into  it.  The  outer  groove  is  used  to 
house  an  0-ring  which  seals  the  fluid  inside  the  chamber.  The  middle  groove  contains  an  0- 
ring  that  seals  the  lower  end  plate  of  the  specimen  to  the  base  plate.  The  inner  groove  seals 
the  lower  cap  drainage  line  to  the  base  plate.  The  lower  end  plate  of  the  soil  specimen  is  3 
in.  in  diameter  and  detachable.  Different  sizes  of  end  plates  can  be  fixed  on  the  base  plate. 

The  apparatus  is  assembled  with  eight  1-1/2  in.  diameter  steel  bolts.  A  specified 
torque,  depending  upon  the  cell  pressure,  is  applied  to  tighten  each  bolt.  Thick  steel  washers 
(0.2  in.)  are  used  under  each  bolt. 


IIL3.4  Pressure  Sources 


Description  of  the  components  used  to  generate  the  high  hydraulic  pressures. 

III.3.5  Pressure  Source  for  Vertical  Loading 

The  pressure  source  is  generated  by  the  linear  motion  of  a  piston  in  a  hydraulic  actua¬ 
tor.  The  pi.  .on  is  pushed  by  a  20  ton  ball  screw  jack  driven  by  a  step  motor.  A  reduction 
gear,  manufactured  by  Sumitomo  SM-Cyclo  model  H3090,  has  a  1:43  gear  ratio.  It  is  con¬ 
nected  between  the  step  motor  and  jack.  The  ball  screw  jack  is  manufactured  by  "Action 
Jack".  The  worm  gear  driven  ball  screw  jack  is  highly  efficient.  It  uses  1/3  the  input  torque 
of  a  machine  screw  jack.  As  a  result  less  horsepower  is  required.  The  ball  screw  jack  used 
has  a  24  in.  travel.  This  computer  controlled  system  can  generate  pressures  up  to  3000  psi. 

UI.3.6  Cell  Pressure  Sources 

Two  pressures  generating  systems  are  used.  One  is  computer  controlled,  the  other  is 
manual. 

The  first  system  is  very  similar  to  the  one  described  above.  The  Sumitomo  gear, 
model  H3090,  has  a  1:17  ratio.  It  is  coupled  with  a  10  ton  ball  screw  action  jack.  The  ball 
screw  jack  has  a  10  in.  travel  distance  and  pushes  an  "Enerpack  RD  910"  hydraulic  cylinder. 
A  controllable  cell  pressure  of  10,000  psi  is  generated  using  this  system. 

The  second  system  uses  an  air  operated  booster  pump.  Manufactured  by  "SC  Pumps 
model  SC40-500-10",  constant  pressures  up  to  20,000  psi  can  be  generated  using  100  psi 
input  air  pressure.  The  pressure  developed  by  the  SC  booster  pump  is  delivered  to  a  "Tescom 
series  54-2000"  hydraulic  pressure  regulator  with  venting  capabilities.  The  pressure  regulator 
has  a  maximum  working  pressure  of  15,000  psi.  It  is  manually  adjusted  using  the  160  MPa 
"Heise"  Bourdon  tube  dial  gauge. 


111,3.7  Control  Panel  and  Pressure  Lines 

All  the  valves  are  made  by  "Snotrik"  and  have  a  capacity  of  45,000  psi.  Two  relief 
valves  are  used:  one,  made  by  "Snotrik"  has  a  range  from  0  to  6,000  psi,  the  other  made  by 
"HIP"  has  range  from  0  to  20,000  psi.  All  high  pressure  fittings  (0  to  20,000  psi)  were 
manufactured  by  Snotrik,  the  low  pressure  fittings  (0  to  3,000  psi)  are  made  by  Swagelok. 

Two  high  precision  pressure  gauges  are  used  to  measure  the  hydraulic  pressure.  A 
"Heise  model  CM-12"  with  a  range  of  0  to  160  MPa  measures  the  chamber  pressure.  It  has 
an  accuracy  of  0.2  MPa.  The  other  gauge  "Heise  model  CMM-12"  has  a  range  of  0  to  40 
MPa,  and  it  is  used  to  measure  the  vertical  pressure  with  an  accuracy  of  0.025  MPa. 

A  high  pressure  flexible  hose  is  used  in  pressurizing  the  cell.  It  has  a  working  pres¬ 
sure  of  25,000  psi.  The  pressure  lines  are  made  from  stainless  steel.  They  are  1/4  in.  O.D. 
with  a  0.065  in.  thickness. 

HI.4  Computerized  Datalogging  and  Test  Control  Systems 
m.4.1  Typical  Process  Flow  for  Controlling  a  Test 

A  typical  test,  whether  drained  or  undrained,  or  whether  under  strain  control  or  stress 
control,  has  very  similar  process  flow  to  control  the  test.  Below  is  described  in  a  simplified 
way  the  interfacing  between  instrumentation,  signal  conditioning,  A/D  conversion,  processing 
raw  data,  the  control  algorithm,  stepping  motor  control,  and  the  hydraulic  syste.a. 

A.  Instrumentation 

The  instrumentation  consists  of  the  load  cell,  which  measures  the  vertical  load  on  the 
triaxial  cell  during  any  type  of  test,  the  LVDT,  which  measures  the  vertical  deflection  of  the 
specimen  during  any  type  of  test,  the  volume  change  device,  which  measures  the  amount  of 
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water  squeezed  from  the  specimen  under  load  during  a  drained  test,  the  pressure  transducers, 
which  measure  the  cell  pressure  during  any  type  of  test,  the  vertical  load  hydraulic  cylinder 
pressure  during  any  type  of  test,  and  the  specimen  pore  pressure  under  load  during  an 
undrained  test.  These  are  all  electrical  devices  whose  output  voltage  changes  as  either  the 
load,  displacement  or  pressure  changes.  The  output  from  these  devices  are  connected  to  the 
signal  conditioning  through  shielded  data  cabling.  All  of  the  instrumentation  is  physically 
located  on  or  near  the  triaxial  cell  on  the  loading  frame. 

B.  Signal  Conditioning 

The  signal  conditioning  performs  several  functions.  It  provides  the  excitation  voltage 
to  drive  the  instrumentation.  The  excitation  voltage  is  10  VDC  for  the  load  cells  and  pressure 
transducers,  and  it  is  about  5.1  VAC  at  about  6.3  kilohertz  for  the  LVDT  and  volume  change 
device.  The  signal  conditioning  also  applies  a  gain  to  the  output  voltage  from  the  instrumenta¬ 
tion  so  it  can  more  fully  occupy  the  input  %  oltage  range  of  the  A/D  converter.  This  maxim¬ 
izes  the  resolution  of  the  A/D  converter.  The  output  of  the  load  cell  and  pressure  transducers 
are  a  maximum  of  20  r  V,  so  a  large  gain  must  be  applied  to  occupy  the  input  voltage  range, 
+-10  VDC,  of  the  A/D  converter.  The  signal  conditioning  system  also  filters  the  incoming 
signal  to  reduce  the  random  higher  frequency  noise  level.  The  output  of  the  signal  condition¬ 
ing  system  is  connected  with  ribbon  cables  to  the  A/D  converter  in  the  microcomputer.  The 
signal  conditioning  is  physically  located  in  a  separate  chassis  near  the  microcomputer. 

C.  A/D  Converter 

The  A/D  converter  converts  the  analog  DC  voltage  signals  to  a  digital  format  that  can 
be  used  by  the  microcomputer.  It  scans  the  six  different  channels  used  by  the  instrumentation 
and  is  connected  to  them  through  the  signal  conditioning.  A  multiplexer  changes  1  .e  channels 
of  the  A/D  converter  and  then  measures  the  voltages.  The  dat .  acquisition  software  that  con- 
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trols  the  A/D  converter,  which  is  in  the  control  program,  then  stores  these  values  in  variable 
numerical  arrays  in  the  microcomputer’s  RAM  memory.  These  values  are  accessible  by  the 
control  software  to  use  for  control  of  the  test. 

D,  Processing  Data  and  Control  Algorithm 

The  software  converts  the  raw  data,  which  are  voltages  from  the  instrumentation,  to 
real  values  by  subtracting  the  reference  voltage  and  multiplying  by  the  calibration  factor  as 
determined  from  the  instrument  calibration. 

Under  strain  control  the  theoretical  strain  or  vertical  displacement  of  the  specimen  is 
known  throughout  the  test.  This  theoretical  value  is  compared  to  the  real  value  measured  by 
the  instrumentation.  The  control  algorithm  computes  a  new  speed  for  the  stepping  motor  that 
controls  the  vertical  movement  to  keep  as  close  to  the  theoretical  value  of  displacement  as 
possible.  Under  stress  control  the  theoretical  vertical  stress  and  theoretical  cell  pressure  are 
known  throughout  the  test.  These  theoretical  values  are  compared  with  the  real  values  as 
measured  by  the  instrumentation.  The  control  algorithm  computes  new  speeds  for  the  step¬ 
ping  motors  that  control  the  vertical  stress  and  the  cell  pressure  to  keep  as  close  to  the 
theoretical  values  as  possible.  The  control  algorithm  computes  the  new  stepping  motor 
speeds.  For  a  more  detailed  description  of  the  software,  see  in. 8  SOFTWARE  DESCRIP¬ 
TION. 

E.  Stepping  Motor  Control 

The  control  algorithm  computes  new  speeds  for  the  stepping  motors.  To  communicate 
this  new  speed,  which  is  in  a  higher-level  computer  language,  to  the  stepping  motors,  a  series 
of  steps  must  be  involved.  First,  if  the  speed  chat  ;e  of  the  stepping  motor  between  the 
current  speed  and  the  prior  speed  is  large  enough  or  if  there  is  a  directional  change,  the  motor 
may  have  to  be  ramped  up  or  down  (make  the  total  speed  change  in  discrete  increments). 
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This  is  performed  in  the  higher-level  computer  language.  The  actual  hardware  calls  are  per¬ 
formed  here.  Second,  the  hardware  calls  are  performed  through  an  assembly  language  inter¬ 
face,  which  ensures  that  the  calls  are  compatible  with  the  higher  -level  language  calling  con¬ 
ventions  and  the  peripheral  interface  device.  Third,  the  above  hardware  calls  access  an 
assembly  language  object  file  with  subroutines  that  control  the  stepping  motor  driver  card 
within  the  microcomputer.  The  stepping  motor  driver  card  sends  out  digital  signals  to  the 
stepping  motor  driver  pack  that  control  stepping  motor  direction,  speed  and  number  of  steps 
to  move.  The  stepping  motor  driver  pack  interprets  these  signals  and  provides  DC  power  to 
the  stepping  motors  in  the  proper  polarities  and  amounts.  For  a  more  detailed  description  of 
the  software,  see  ffl.8  SOFTWARE  DESCRIPTION. 

F.  Hydraulic  System 

Each  stepping  motor  is  connected  to  the  hydraulic  system  through  a  reduction  gear. 
The  reduction  gear  drives  a  ball-screw  jack,  which  also  has  some  reduction  effect  The  action 
jack  converts  rotational  motion  into  linear  motion.  The  action  jack  pushes  or  pulls  on  a 
hydraulic  cylinder,  which  injects  or  withdraws  hydraulic  fluid  into  the  hydraulic  system. 

One  stepping  motor/hydraulic  system  injects  or  withdraws  hydraulic  fluid  from  the 
loading  frame’s  vertical  load  cylinder.  The  vertical  load  cylinder  is  located  in  the  base  of  the 
loading  frame  and  actuates  a  piston  that  moves  in  the  vertical  direction.  This  hydraulic 
cylinder  piston  has  a  small  platform  on  which  the  triaxial  cell  rests.  The  triaxial  cell  piston  is 
pressed  against  the  loading  frame’s  stationary  crosshead  and  moves  relative  to  the  hydraulic 
cylinder  piston.  Therefore,  this  stepping  motor/hydraulic  system  controls  the  vertical  move¬ 
ment  or  vertical  load  of  the  specimen  in  the  triaxial  cell. 

The  other  stepping  motor/hydraulic  system  injects  or  withdraws  hydraulic  fluid  from 
the  triaxial  cell  to  directly  increase  or  decrease  the  cell  pressure  on  the  specimen  in  the  triax- 
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ial  cell.  The  system  instrumentation  is  connected  to  the  triaxial  cell,  so  the  process  flow 
would  go  back  to  A.  INSTRUMENTATION  and  loop  indefinitely  until  the  end  of  the  test. 

The  above  description  of  the  hydraulic  system  does  not  describe  all  of  the  piping, 
shutoff  and  isolation  valving,  and  gages. 

m.5  System  Performance  and  Control  Capabilities 

The  software  limits  the  speed  of  the  stepping  motors  to  a  maximum  value  of  2,500 
steps  per  second  in  half-step  mode.  This  is  equivalent  to  6.25  revolutions  per  second.  The 
limit  is  intended  to  prevent  motor  stalling  under  maximum  loads.  This  maximum  speed 
translates  into  a  maximum  deformation  rate  of  about  3  inches  per  hour.  The  maximum  stress 
rate  for  the  cell  pressure  stepping  motor/hydraulic  system  is  7.46  MPa  per  minute,  assuming 
that  the  motor  speed  is  2,500  steps  per  second  and  that  the  specimen  has  the  same  stiffness  as 
the  system.  The  maximum  stress  rate  for  the  deviator  stress/hydraulic  system  is  2.35  MPa  per 
minute,  assuming  that  the  motor  speed  is  2,500  steps  per  second  and  that  the  specimen  has  the 
same  stiffness  as  the  system. 

The  random  noise  level  varies  depending  on  the  instrumentation  involved  and  the  gain 
applied  to  the  output  voltage  from  the  instrumentation.  They  are: 

1.  Load  Cell  -  0.008 %  error.  The  gain  is  the  same  for  the  20  ton  and  the  50-ton  load 
cells.  Therefore,  the  0.008%  error  represents  1.6  kilogram  for  the  20-ton  load  cell  and 
4.0  kilograms  for  the  50-ton  load  cell. 

2.  Pressure  Transducers  -  0.06%  error.  The  gain  will  vary  depending  on  the  pressure 
range  selected.  The  current  range  is  at  8  MPa.  Therefore,  this  error  represents  a  value 
of  0.005  MPa  or  0.7  pounds  per  square  inch. 
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3.  Volume  Change  Device  -  0.003%  error.  The  gain  will  not  change  on  this  device. 

Therefore,  its  error  is  equivalent  to  0.003  cubic  centimeters. 

4.  LVDT  -  0.004%  error.  The  gain  will  not  change  on  this  device.  Therefore,  its  error  is 

equivalent  to  0.0001  inch. 

The  above  values  for  random  noise  in  the  system  do  not  reflect  the  true  total  error  in 
the  instrumentation,  because  each  piece  of  instrumentation  has  its  own  accuracy.  See  HI. 7 
EQUIPMENT  DESCRIPTION  for  hardware  instrumentation  tolerances.  In  most  cases  the 
accuracy  of  the  system  will  be  the  hardware  accuracy,  since  the  random  noise  level  is  much 
less  than  the  hardware  tolerance. 

The  control  program’s  ability  to  control  the  hydraulic  pressures  and  deformations  in 
the  system  are  based  more  on  the  random  noise  level  than  the  overall  accuracy  of  the  system, 
since  the  random  noise  represents  the  repeatability  of  the  data.  The  overall  system  accuracy 
is  based  more  on  the  hardware  tolerance,  since  it  represents  the  deviations  of  measured  data 
and  reality.  Therefore,  the  following  values  representing  the  ability  for  controlling  parameters 
are  based  upon  the  random  noise  levels,  rather  than  the  hardware  tolerance.  The  overall  sys¬ 
tem  accuracy  is  still  the  hardware  tolerance.  This  is  emphasized  to  allow  the  control  software 
to  be  evaluated  independently  of  the  instrumentation  hardware.  This  is  important,  because  if 
the  current  instrumentation  is  replaced  by  higher-accuracy  instrumentation,  then  th-3  control 
software  has  the  ability  to  improve  the  overall  system  accuracy  to  a  level  of  the  values  shown 
below. 


The  control  algorithm  for  strain  control  with  constant  cell  pressure  can  control  the 
strain  rate  to  within  0.1%  of  the  strain  rate  and  control  the  cell  pressure  to  a  level  of  about 
0.2%  of  the  cell  pressure  when  it  is  at  full  scale  voltage,  10  VDC.  The  control  algorithm  for 
stress  control  can  control  the  cell  pressure  and  deviator  stress  to  about  the  same  level,  but 
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there  is  an  additional  factor  to  consider:  the  triaxial  cell  piston  seal  friction.  The  triaxial  cell 
piston  seal  friction  appears  to  be  a  function  of  cell  pressure.  For  the  standard  triaxial  tests 
under  constant  cell  pressure,  the  dynamic  friction  can  be  easily  determined  by  the  software, 
since  its  direction  of  motion  is  only  in  one  direction  and  the  test  is  performed  at  constant  cell 
pressure.  However,  for  a  constantly  varying  cell  pressure  with  changes  of  direction  in  vertical 
motion,  it  is  more  difficult.  The  model  used  in  the  control  software  can  add  an  additional 
maximum  error  of  approximately  0.4%,  based  on  the  maximum  deviator  stress  for  the  type  of 
soils  currently  used.  This  error  occurs  at  vertical  motion  changes  in  direction.  Points  that  are 
not  located  near  vertical  motion  changes  in  direction  appear  to  have  a  much  lower  additional 
error  of  about  0.1%  to  0.2%.  The  selection  of  points  in  the  stress  path  is  crucial  to  keep  this 
additional  error  to  a  minimum.  If  the  user  can  keep  the  areas  of  the  stress  path  he  wishes  to 
examine  away  from  the  vertical  motion  changes  in  direction,  he  can  generally  keep  this  addi¬ 
tional  error  to  a  level  of  about  0. 1  %. 

III.6  System  Safety 

Obviously,  operating  under  such  high  pressures  presents  significant  safety  hazards  to 
the  operators  in  the  event  of  equipment  failure.  Therefore  system  safety  is  designed  with 
several  levels  of  protection,  both  in  hardware  and  software. 

III.6.1  Hardware  Safety  Features 

The  hardware  is  comprised  of  the  following: 

1.  Limit  Switches  -  Both  hydraulic  cylinders  driven  by  the  stepping  motors  have  limit 
switches  installed  to  limit  the  maximum  length  of  travel.  The  limit  switches  switch  off 
the  power  to  the  stepping  motor  driver  pack,  which  stops  the  stepping  motors. 

2.  Pressure  Relief  Valves  -  Both  high  pressure  stepping  motor/hydraulic  systems  have 
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pressure  relief  valves  installed  to  limit  the  maximum  pressure  in  each  system.  They  are 
preset  to  16,500  pounds  per  square  inch  for  the  cell  pressure  system  and  3,600 
pounds  per  square  inch  for  the  vertical  load  system.  There  is  also  a  pressure  relief 
valve  on  the  volume  change  device/backpressure  system  preset  to  100  pounds  per 
square  inch. 

3.  Physical  Shield  -  There  is  a  one-inch  thick  steel  plate  measuring  18  inches  by  47 
inches  mounted  on  the  loading  frame  to  help  protect  the  operator,  if  the  triaxial  cell 
fails. 

4.  Monitoring  Devices  -  There  are  direct  manual  monitoring  devices  to  allow  the  operator 
to  monitor  system  parameters  during  the  test.  They  consist  of  precision  pressure  gages 
and  dial  gages.  If  the  parameters  reach  levels  beyond  the  system  ratings,  the  operator 
can  manually  stop  the  test. 

5.  Emergency  Stop  Switch  -  There  is  an  emergency  stop  switch  that  directly  issues  an 
uninterruptible  stop  order  to  the  stepping  motor  driver  card,  which  controls  the  step¬ 
ping  motors.  The  physical  location  of  the  emergency  stop  switch  is  on  the  signal  con¬ 
ditioning  chassis  located  next  to  the  microcomputer. 

III.6.2  Software  Safety  Features 

1.  Control  Program  -  The  control  program  monitors  the  system  instrumentation  and  has 
software  routines  checking  the  vertical  load,  the  cell  pressure,  and  the  loading  frame 
vertical  load  hydraulic  cylinder  pressure  to  ensure  that  they  remain  withing  the  system 
ratings.  If  not,  the  software  will  issue  a  visual  and  audible  warning  and  terminate  the 
test. 
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The  test  software  constantly  monitors  the  test  progress  and  displays  four  plots  on  the 

screen  in  real  time  to  aid  the  operator  in  monitoring  the  test  progress.  During  the  test,  the 

operator  can  stop  the  test  at  any  time  by  pressing  a  special  function  key  on  the  microcomputer 

keyboard. 

III.7  Equipment  Description 

A.  Instrumentation 

1.  Pressure  Transducers  -  HBM,  Model  P  3MA,  1000  bar  and  500  bar  capacities,  2 
mV/V  sensitivity,  0.15%  tolerance. 

2.  Load  Cells  -  HBM,  Model  C3H2,  20  metric  ton  and  50  metric  ton  capacities,  2  mV/V 
sensitivity,  0.1%  tolerance. 

3.  LVDT  -  Schaevitz  Engineering,  Model  2000  E,  linearity  0.56%,  2  inch. 

4.  Volume  Change  Measuring  Device  -  Custom  made  device  using  a  double-acting  piston, 
which  alternately  pushes  or  expells  water  from  one  of  two  chambers  through  a  four¬ 
way  valve.  An  LVDT  is  attached  to  the  piston  to  measure  movement. 

5.  Cabling  -  Shielded,  4  -  wire,  stranded  type. 

B.  Signal  Conditioning 

1.  System  -  Analog  Devices,  3B  Series  Signal  Conditioning  Subsystem,  AC/DC  power 
supply,  terminal  connections  for  8  separate  channels  of  inputs. 

2.  Signal  Conditioning  Modules  -  Analog  Devices,  Model  3B16  for  strain  gage  input, 
adjustable  gain,  adjustable  zero,  10  VDC  excitation  voltage,  module  modified  to  allow 
calibration  grade  potentiometer  to  vary  gain.  Model  3B17  for  LVDT  input,  adjustable 
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zero,  adjustable  span,  adjustable  gain,  adjustable  AC  excitation  voltage,  adjustable 
excitation  frequency. 

3.  Signal  Conditioner! AD  Converter  Interface  Board  -  Prototype  board  to  interface  the 
signal  conditioning  system,  which  is  a  single-ended  input  type  to  the  AD  converter, 
which  is  a  differential-input  type. 

C.  A/D  Converter 

Data  Translation,  Model  DT2801/5716,  16  bit  resolution,  2.5  kilohertz  maximum 
conversion  speed,  8  input  channel  capacity,  differential  input,  +-10  VDC  maximum  voltage 
input. 

D.  Microcomputer 

IBM,  Model  AT,  512  kilobytes  of  RAM,  6  megahertz  clock  speed,  1.2  megabyte  flexi¬ 
ble  disk  drive,  20  megabyte  fixed  disk  drive. 

E.  Stepping  Motor  Driver  Card 

Rogers  Labs,  Model  I  432  I/O,  4-axis  motion  stepping  motor  capability. 

F.  Stepping  Motor  Interface  Board 

Termination  board  between  stepping  motor  driver  card  and  stepping  motor  driver  pack. 

G.  External  Relay  Board 

Relay  board  which  allows  external  devices  such  as  switches  to  provide  direct  input  to 
the  stepping  motor  driver  card.  Three  DC  relays  and  AC/DC  power  supply. 
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H.  Stepping  Motors 


Superior  Electric,  Model  M111-FD12,  4  phase,  2.26  VDC,  6.1  amps  per  phase. 
Anaheim  Automation,  Model  42D112-12AR,  4  phase,  3.6  VDC,  6.1  amps  per  phase. 

I.  Stepping  Motor  Driver  Pack 

Anaheim  Automation,  Model  DPB05,  4  phase,  6.5  amps  per  phase  capacity,  bilevel 
type,  two  motor  output. 

J.  Power  Conditioner 

Topaz,  Model  Peak  2,  02406-01P3,  120  VAC,  1000  VA  capacity. 

K.  Surge  Protector 

Inmac,  Model  8207,  4  outlet,  120  VAC,  15  amp  rating,  zener  diode,  MOV,  gas-tube 
protection  levels. 

ni.8  Software  Description 

There  are  several  programs  that  were  written  for  different  purposes.  All  the  higher- 
level  language  programs  were  written  using  Microsoft  QuickbasicR  Version  3.0.  All  low 
level  language  programs  were  written  in  Assembly  Language  and  assembled  using  the  Micro- 
softR  Macro  Assembler  Version  5.0. 

The  programs  fall  into  two  general  groups,  utility  programs  and  testing  programs. 

A.  Utility  Programs 

1.  MENU.EXE  -  This  program  provides  a  user  interface  for  inputting  all  of  the  data 
relevant  to  run  the  desired  test  It  is  fully  menu  driven  and  prompts  the  user  for  infor- 


mation  and  also  has  error  checking  routines.  The  test  data  is  written  into  an  ASCII 
file  for  access  by  the  test  programs.  The  program  length  is  approximately  1,300  lines. 

2.  DATAOUT.EXE  -  This  program  provides  the  user  with  the  amplified  instrumentation 
voltages  on  the  different  data  channels,  whenever  it  is  triggered.  The  voltages  are 
written  to  an  ASCII  file.  This  is  used  when  calibrating  the  instrumentation.  The  pro¬ 
gram  length  is  approximately  200  lines. 

3.  STEPPER.EXE  -  This  program  allows  the  user  to  independently  operate  the  two  step¬ 
ping  motors  in  either  direction  and  at  any  speed.  The  program  length  is  approximately 
200  lines. 

4.  QT5DATA.EXE  -  This  program  provides  the  user  with  the  amplified  instrumentation 
voltages  on  the  different  data  channels,  whenever  it  is  triggered  and  is  similar  to 
DATAOUT.EXE,  except  it  is  utilized  with  the  stepping  motor  operating.  This  allows 
the  user  to  calibrate  in  a  continuous  running  mode,  so  hysteretic  behavior  of  instru¬ 
mentation  can  be  checked.  Voltages  are  written  to  an  ASCII  file.  The  program  length 
is  approximately  400  lines. 

5.  FRICT.EXE  -  This  program  allows  the  user  to  measure  the  triaxial  cell  piston  seal 
static  and  dynamic  friction  forces.  Its  method  of  operation  is  to  increase  the  triaxial 
cell  pressure  to  a  series  of  desired  values  and  then  run  the  piston  both  up  and  down, 
while  keeping  the  cell  pressure  constant.  The  friction  force  is  the  difference  between 
the  converted  values  of  the  load  cell  and  the  cell  pressure  multiplied  by  the  piston 
area.  Therefore,  this  program  measures  both  the  static  friction  (when  the  piston  first 
starts  to  move)  and  dynamic  friction(after  the  piston  has  moved  several  thousandths  of 
an  inch)  force.  All  values  are  written  to  an  ASCII  file.  The  program  length  is  approx¬ 
imately  700  lines. 
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6.  FR1CT1.EXE  -  This  program  allows  the  user  to  measure  the  piston  seal  dynamic  fric¬ 
tion  force.  This  program  is  different  from  FRICT.EXE,  because  it  does  not  stop  at  a 
series  of  cell  pressures.  It  continuously  increases  the  cell  pressures  up  and  down  to 
desired  values  and  takes  continuous  readings  of  the  friction  force.  Since  the  cell  pres¬ 
sure  is  increasing  continuously,  the  value  that  is  measured  is  the  dynamic  friction 
force.  All  values  are  written  to  an  ASCII  file.  The  program  length  is  approximately 
500  lines. 

7.  FRICT2.EXE  -  This  program  is  similar  to  FRICT1.EXE  except  it  has  incorporated 
within  it  the  subroutine  that  predicts  triaxial  cell  seal  friction,  so  the  real  friction  can 
be  compared  co  the  estimated  friction  from  the  friction  model  used  in  the  test  pro¬ 
grams.  The  program  length  is  approximately  600  lines. 

B.  Test  Programs 

1.  TEST.EXE  -  This  program  is  the  access  program  for  all  the  different  test  programs. 
The  user  executes  this  program  and  it  checks  the  data  file  created  by  the  menu  pro¬ 
gram  and  automatically  executes  the  correct  test  program.  The  program  length  is 
approximately  20  lines. 

2.  CU_CD_S3.EXE  -  This  program  controls  two  different  types  of  tests.  They  are  the 
drained  and  undrained  triaxial  compression  tests  with  constant  cell  pressure  performed 
under  constant  strain  rate.  Isotropic  and  anisotropic  consolidation  are  also  included. 
See  below  for  more  detailed  description  of  program  elements.  The  program  length  is 
approximately  2,400  lines. 

3.  CU_CD_S1  .EXE  -  This  program  controls  two  different  types  of  tests.  They  are  the 
drained  and  undrained  triaxial  compression  tests  with  constant  vertical  stress  performed 
under  constant  strain  rate.  Isotropic  and  anisotropic  consolidation  are  included.  See 
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below  for  more  detailed  description  of  program  elements.  The  program  length  is 
approximately  1,800  li  TICS. 

4.  GENSTRSS.EXE  -  This  program  controls  a  general  stress  path  type  of  test  performed 
under  stress  control.  See  below  for  more  detailed  description  of  program  elements. 
The  program  length  is  approximately  1,900  lines. 

5.  GENSTRN.EXE  -  This  program  controls  a  general  strain  path  type  of  test  performed 
under  strain  control.  See  below  for  more  detailed  description  of  program  elements. 
The  program  length  is  approximately  2,100  lines. 

C.  Test  Program  Elements 

The  different  test  programs  are  composed  of  several  different  elements.  The  major 
elements  are  data  acquisition,  data  reduction,  storing  data,  control,  actuation  of  stepping 
motors,  and  graphic  display  of  test  data. 

1.  DATA  ACQUISITION  ELEMENT  -  The  data  acquisition  element  is  written  in  the  basic 
computer  language,  since  the  A/D  converter  is  designed  to  accept  commands  using  the 
"memory  mapped  I/O"  concept,  which  uses  different  RAM  memory  addresses  as  Com¬ 
mand,  Data  and  Error  registers.  The  basic  computer  language  is  adept  at  this  function. 
Commands  are  sent  to  the  Command  register’s  memory  address  in  specified  bit  pat¬ 
terns.  Data  is  read  from  the  Data  register  address  in  two-byte  increments  and  con¬ 
verted  from  binary  to  decimal.  Error  messages  are  read  from  the  Error  register  address 
and  their  meaning  is  translated  by  their  bit  patterns.  This  program  element  samples 
data  at  a  rate  of  approximately  2  kilohertz.  Six  channels  are  monitored,  which 
correspond  to  the  vertical  load  cell,  the  specimen  pore  pressure,  the  triaxial  cell  pres¬ 
sure,  the  loading  frame’s  vertical  load  hydraulic  cylinder  pressure,  the  volume  change, 
and  the  vertical  motion  LVDT.  The  data  received  from  the  A/D  converter  is  a  DC 
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signal  voltage  from  the  different  pieces  of  instrumentation.  Fifty  samples  on  each 
channel  are  taken  and  then  averaged  to  act  as  a  low-pass  data  filter,  which  tends  to 
reduce  the  random  noise  level. 

2.  STORING  DATA  -  The  above  data  is  stored  onto  the  fixed  disk  in  an  ASCII  file  called 
VOLTAGE.PRN.  The  individual  channel  voltages  are  then  converted  into  their  real 
instrumentation  values  by  subtracting  the  reference  voltage  and  multiplying  by  the  cali¬ 
bration  factor  for  that  piece  of  instrumentation.  These  "real"  values  are  then  stored 
onto  fixed  disk  in  an  ASCII  file  called  DATA.PRN.  When  the  actual  shearing  portion 
of  the  test  starts,  these  values  from  the  instrumentation  are  converted  to  actual 
engineering  parameters  for  axial  strain,  volumetric  strain,  deviator  stress,  stress  ratios, 
and  P  and  Q,  which  are  stress  path  parameters.  Control  is  accomplished  using  these 
values.  These  engineering  parameters  are  then  stored  onto  fixed  disk  in  an  ASCII  file 
called  FINAL.PRN.  Since  the  three  above  listed  data  files  are  ASCII  files,  they  can  be 
imported  directly  into  LOTUS  1-2-3R  for  post-test  plotting  of  results  or  further  data 
analysis. 

If  the  measured  value  from  the  instrumentation  exceeds  specified  values  for 
maximum  load  or  pressure,  the  program  stops  execution  and  automatically  starts 
unloading  the  specimen. 

3.  CONTROL  -  There  are  two  general  types  of  control  elements.  One  for  strain  control 
and  one  for  stress  control. 

Under  strain  control,  at  any  time  during  the  shearing  portion  of  the  test  there 
should  be  a  specific  vertical  displacement  of  the  specimen.  This  can  be  calculated 
from  the  elapsed  time  of  the  test  and  the  strain  rate.  This  theoretical  value  is  com¬ 
pared  with  the  real  value  given  from  the  LVDT  and  if  it  is  less  or  greater  than  this 
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value,  a  new  speed  for  the  vertical  direction  stepping  motor  is  computed.  The  control 
algorithm  uses  the  method  of  proportional  response.  In  other  words  it  adjusts  the 
change  in  stepping  motor  speed  relative  to  the  amount  of  continuing  divergence  from 
the  theoretical  displacement  value.  It  also  buffers  this  proportional  response  so  that  if 
the  difference  between  the  theoretical  and  real  values  of  displacement  is  closing  at  an 
acceptable  rate  it  does  not  change  the  stepping  motor  speed.  The  effect  of  this  is  to 
eliminate  overshooting  the  theoretical  value  and  help  keep  the  control  algorithm  stable. 

Some  of  the  tests  require  constant  cell  pressure.  Since  the  piston  is  descending 
into  the  triaxial  cell,  fluid  needs  to  be  removed  from  the  triaxial  cell  to  keep  a  constant 
cell  pressure.  The  other  stepping  motor  bleeds  off  this  excess  fluid,  however,  it  is  con¬ 
trolled  by  using  the  measured  value  of  cell  pressure,  comparing  it  to  the  desired  cell 
pressure.  This  control  algorithm  is  very  similar  to  the  above  one  for  strain  control. 

Under  stress  control,  at  any  time  during  the  test  there  should  be  a  specific  verti¬ 
cal  stress  and  a  specific  cell  pressure  on  the  specimen.  This  can  be  calculated  from 
the  elapsed  time  of  the  test  and  the  vertical  and  horizontal  stress  rate.  These  theoreti¬ 
cal  values  are  compared  with  the  real  values  given  from  the  cell  pressure  transducer 
and  a  calculated  value  involving  the  load  cell,  cell  pressure,  seal  friction,  and  mem¬ 
brane  stiffness.  If  the  real  values  are  less  or  greater  than  these  theoretical  values,  new 
speeds  for  the  vertical  direction  stepping  motor  and  the  cell  pressure  stepping  motor 
are  computed.  The  control  algorithm  also  uses  the  method  of  proportional  response. 
In  other  words  it  adjusts  the  change  in  stepping  motor  speed  relative  to  the  amount  of 
continuing  divergence  from  the  theoretical  vertical  stress  and  cell  pressure  values.  It 
also  buffers  this  proportional  response  so  that  if  the  difference  between  the  theoretical 
and  real  values  of  vertical  stress  and  cell  pressure  are  closing  at  an  acceptable  rate,  it 
does  not  change  the  stepping  motor  speed.  The  effect  of  this  is  to  eliminate 
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overshooting  the  theoretical  values  and  helps  keep  the  control  algorithm  stable. 

4.  ACTUATION  OF  STEPPING  MOTORS  -  This  program  element  operates  the  stepping 
motors.  It  is  composed  of  three  parts.  They  are  the  higher-level  language  part,  the 
assembly  language  interface,  and  the  vendor-supplied  proprietary  object  file. 

The  vendor-supplied  proprietary  object  file  is  a  series  of  assembled  assembly 
language  subroutines,  which  can  be  called  from  higher-level  languages  through  an 
interface.  Some  of  the  different  subroutines  are  FREQ  1N2 -and  FREQ3N4,  which  set 
the  desired  speed  of  up  to  four  stepping  motors,  MOT1,  MOT2,  MOT3,  and  MOT4, 
which  set  the  direction  of  rotation  of  the  stepping  motors  and  the  number  of  steps  to 
travel  in  that  direction.  There  are  some  limits  to  the  control  of  the  stepping  motors. 
Of  particular  importance  is  the  limitation  that  the  stepping  motor  speed  must  be  above 
40  steps  per  second  and  that  a  motor  can  only  be  commanded  to  move  between  1  and 
32767  steps  in  one  subroutine  call.  The  significance  of  the  40  steps  per  second  limita¬ 
tion  is  that  if  the  constant  stepping  motor  speed  is  directed  to  fall  below  40  steps  per 
second  by  the  control  algorithm,  then  the  only  way  to  approximate  this  low  speed  is  to 
operate  the  motor  in  an  intermittent  (on/off)  mode  instead  of  at  a  continuous  speed. 
The  significance  of  the  1  to  32767  steps  limitation  is  that  when  it  is  desired  to  run  the 
motor  continuously,  the  motor  must  be  called  repeatedly  before  the  number  of  steps 
from  the  last  subroutine  call  runs  out  and  stalls  the  stepping  motor. 

The  assembly  language  interface,  called  QB1,  was  written  to  provide  a  bridge 
between  the  higher-level  language  and  the  vendor- supplied  assembly  language  subrou¬ 
tines.  This  is  necessary  to  accommodate  the  calling  conventions  of  the  particular 
higher-level  language,  which  is  Microsoft  QuickbasicR  Version  3.0. 
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The  higher-level  language  portion  consists  of  two  elements.  They  are  a  subrou¬ 
tine  called  STEPMOT  and  a  subprogram  called  STEPPER.  They  perform  several 
necessary  functions  that  cannot  be  performed  directly  by  the  vendor-supplied  assembly 
language  subroutines.  STEPMOT  maintains  a  record  of  the  last  and  current  stepping 
motor  directions  and  speeds,  which  is  used  by  STEPPER  to  ramp  up  or  ramp  down 
(make  the  total  speed  change  in  smaller  discrete  steps)  the  stepping  motor  speeds  to 
the  desired  level.  Ramping  is  required  because  stepping  motors  cannot  change  speeds 
in  large  increments  without  stalling.  STEPPER  makes  the  -actual  calls  to  the  vendor- 
supplied  assembly  language  object  file. 

GRAPHIC  DISPLAY  -  This  portion  is  composed  of  two  separate  elements.  They  are 
the  background  screen  generation  element  and  the  line  plotting  element.  Together  they 
display  in  real-time  the  engineering  parameters  calculated  earlier.  Both  elements  were 
written  for  an  Enhanced  Graphics  Adapter  type  of  display.  The  background  element 
splits  the  screen  into  four  separate  graphics  windows.  Each  window  is  for  an  X-Y 
graph  plot  of  deviator  stress  versus  axial  strain,  effective  stress  ratio  versus  axial  strain, 
stress  path  (Cambridge  P  versus  Cambridge  Q),  and  either  volumetric  strain  versus 
axial  strain  or  pore  pressure  versus  axial  strain  depending  on  whether  the  test  is 
drained  or  undrained.  The  vertical  scale  is  auto-scaling  as  a  function  of  the  confining 
pressure.  Window  borders,  axes,  grid  lines  and  labels  are  plotted  on  the  screen.  The 
background  element  is  generated  only  once,  so  when  the  test  is  running  only  lines 
between  data  points  will  need  to  be  plotted  on  the  existing  background. 

The  line  plotting  element  plots  lines  between  the  last  data  point  and  the  current 
data  point  Also,  the  current  numerical  values  of  the  calculated  engineering  parameters 
are  printed  on  the  screen  within  the  window  boundaries. 
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III.9  Results  of  Two  Experiments 

Triaxial  equipment  already  available  in  the  Soil  Mechanics  Laboratory  has  been  used 
to  perform  a  series  of  baseline  tests  in  triaxial  compression  and  extension.  Drained  tests  were 
performed  with  constant  cell  pressures  up  to  600  psi  in  compression  and  extension.  The 
granular  material  used  for  these  tests  and  the  tests  in  the  high  pressure  triaxial  apparatus  is 
Cambria  Sand  of  similar  composition  to  that  used  in  two  previous  research  projects.  This  is 
the  same  sand  that  is  being  modeled  in  the  numerical  investigation  being  conducted  using  the 
PROBS  computer  program. 

The  results  of  two  tests  performed  in  the  high  pressure  triaxial  apparatus  are  shown  in 
the  diagrams  on  the  following  pages.  One  test  is  a  consolidated  -  drained  test  (CD  Test),  and 
the  other  is  a  consolidated  -  undrained  test  (CU  Test).  The  specimens  were  isotropically  con¬ 
solidated  to  1200  psi  (8.277  MPa)  in  both  tests  before  shearing  was  initiated. 

The  results  of  the  drained  test  are  shown  in  Figs.  2,3,4  and  5.  Figure  2  indicates  the 

f 

effective  stress-path  on  a  Cambridge  p^g  diagram.  For  a  drained  test  at  constant  confining 
pressure  the  stress-path  is  a  straight  line  at  a  slope  of  1  horizontal  on  3  vertical.  The  deviator 
stress  -  axial  strain  curve  is  shown  in  Figure  3  and  the  max.  deviator  stress  is  indicated  at 
19.79  MPa  corresponding  to  an  effective  friction  angle  of  33.0°.  The  variation  of  the 
effective  stress  ratio  is  shown  in  Figure  4,  and  the  volume  change  curve  is  shown  in  Figure  5. 
At  this  confining  pressure,  the  initially  dense  soil  specimen  is  highly  compressive,  and  this 
compression  is  directly  related  to  the  crushing  of  the  soil  grains. 

The  results  of  the  undrained  test  are  shown  in  Figures  6,7,8  and  9.  The  effective 
stress-path  in  Figure  6  is  controlled  by  the  variation  in  pore  pressure  as  the  deviator  stress  is 
increased.  The  deviator  stress-axial  strain  variation  shown  in  Figure  7  and  Figure  8  indicates 
the  variation  of  the  effective  stress  ratio.  The  maximum  effective  stress  ratio  of  3.26  is 
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reached  at  about  18%  axial  strain.  This  corresponds  to  a  friction  angle  of  32.0°.  The  pore 
pressure  variation  is  shown  in  Figure  9. 

The  diagrams  of  test  results  shown  in  Figures  2-9  indicate  the  capability  of  the  experi¬ 
mental  equipment  to  perform  tests  of  the  types  required  for  this  investigation. 

IIL10  Microphysical  Characteristics  of  Granular  Materials 

Part  of  the  experimental  study  focuses  on  an  investigation  of  the  microphysical  proper¬ 
ties  of  granular  materials  and  their  changes  during  shear  tests.  This  will  have  important  bear¬ 
ing  on  the  development  of  realistic  computer  models  of  granular  materials,  and  provide 
another  basis  for  comparing  experiments  and  numerical  studies.  In  essence  this  work  involves 
counting  the  number  of  contact  points  and  their  orientation  at  various  points  of  a  shear  test. 
In  particular,  the  evolution  of  microphysical  geometry  of  the  granular  materials  is  being  stu¬ 
died  by  examining  triaxial  specimens  at  several  stages  of  a  test,  including  before,  during  and 
after  failure  of  the  specimen.  Parameters  such  as  shape,  size  distribution,  breakage,  and  orien¬ 
tation  of  particles  as  well  as  the  number  and  orientation  of  particle  contacts  are  being  meas¬ 
ured.  By  "freezing"  the  specimens  at  different  stages,  these  parameters  may  be  measured  and 
their  evolution  can  be  followed  during  the  tests. 

It  was  originally  planned  that  this  part  of  the  project  would  be  performed  on  the  same 
Cambria  Sand  as  used  in  the  other  parts  of  the  project.  However,  it  quickly  became  clear  that 
the  particles  were  too  small  (in  the  range  from  1.0  to  2.0  mm)  for  observation  of  contact 
points  and  their  directions.  This  portion  of  the  study  is  therefore  being  performed  on  another 
size  range  of  Cambria  Sand.  The  grain  sizes  in  this  material  range  from  2.0  mm  to  5.0  mm. 
The  grain  shapes,  distribution,  etc.  in  this  range  have  already  been  studied  and  adjusted 
(where  possible)  to  match  the  original  Cambria  Sand. 
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The  experimental  equipment  for  the  study  of  grain  contact  angles  and  their  evolution 
has  been  constructed  and  refined  for  use  on  a  routine  basis.  Several  specimens  have  been  stu¬ 
died.  Approximately  200  grains  (40  grains  in  each  of  five  layers)  located  near  the  middle  in 
each  of  5-7  specimens  are  sheared  in  triaxial  compression  to  different  amounts  of  vertical 
deformation,  "frozen"  (using  a  paint  technique),  and  taken  apart,  grain  by  grain,  noting  posi¬ 
tion,  orientation,  grain  contact  numbers  and  directions.  By  combining  the  results  from  the  5-7 
specimens  it  is  possible  to  follow  the  evolution  of  these  microphysical  characteristics. 

This  study  is  being  performed  on  specimens  sheared  at  a  low  confining  pressure  of  1.0 
kg/cm  .  Thus,  particle  breakage  is  unlikely  to  occur  in  these  tests.  Particle  breakage,  conse¬ 
quent  variation  in  size  distribution  and  possibly  the  variation  in  shape  (due  to  breakage)  will 
be  studied  in  tests  at  very  high  confining  pressures  as  explained  above. 

ni.ll  Summary  and  Conclusions 

It  took  approximately  the  first  year  of  the  two  year  contract  to  obtain  a  loading 
machine  for  the  high  pressure  triaxial  testing  to  be  performed  on  this  project.  The  reasons  for 
this  time  delay  have  been  explained  in  the  progress  reports  and  in  the  present  report.  How¬ 
ever,  at  this  time,  i.e.  at  the  end  of  the  two-year  period,  the  experiment  facility  for  high  pres¬ 
sure  triaxial  testing  has  been  constructed  and  it  is  in  excellent  working  condition.  The  capa¬ 
bilities  of  the  equipment  meets  or  exceeds  the  originally  intended  levels  of  performance,  and 
it  has  been  constructed  within  the  original  budget.  The  major  cost  has  been  the  extra  time 
spent  in  reaching  the  point  where  the  loading  frame  was  in  house  at  UCLA. 

High  pressure  triaxial  testing  is  progressing  according  to  a  regular  schedule  set  up  to 
provide  data  to  be  analyzed  in  relation  to  stability/instability  of  granular  soils  at  high  pres¬ 
sures.  It  is  anticipated  that  the  results  of  these  tests  and  their  analyses  will  be  presented  in 
addendums  to  the  present  final  report.  Along  with  these  results  will  also  be  more  detailed 
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presentations  of  the  experimental  equipment  as  well  as  the  study  of  micromechanical  charac¬ 
teristics  of  granular  materials  initially  proposed  to  AFOSR. 
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200  kips  Gilmore  Loading  Machine 
(4  columns) 


HI- 1-  Schematic  Diagram  of  Experimental  Setup  for  High  Pres¬ 
sure  Triaxial  Testing. 


Triaxia  1  Compression  of  Cambria  Sand 

CD  Test;  Eff.  Conf.  Pressure-  1200  psi 
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III.2  Effective  Stress-Path  for  Drained  Test  on  Cambria  Sand. 


Triaxial  Compression  of  Cambria  Sand 

CD  Test;  Eff.  Conf.  Pressure*3  1200  psi 


III. 3  Deviator  Stress-Axial  Strain  Relation  for  Drained  Test  on 
Cambria  Sand. 
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Triaxial  Compression  of  Cambria  Sand 

CD  Test;  Eff.  Conf.  Preasure=  1200  psi 
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Triaxial  Compression  of  Cambria  Sand 

CU  Test;  Eff.  Conf.  Pressure  -•  1200  psi 


III.6  Effective  Stress-Path  in  p!-q  Undrained  Test  on  Cambria 
Sand. 


Triaxia I  Compression  of  Cambria  Sand 

CU  Test;  Eff.  Conf.  Pressure  »  1200  psi 


Triaxial  Compression  of  Cambria  Sand 

CU  Test;  Eff.  Conf.  Pressure  =  1200  psi 
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Pressure  Variation  in  Undrained  Test  on  Cambria 
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